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An  extensive  analysis  that  explicitly  illustrates  how  speckle  effects 
interfere  with  COAT  operation  has  been  developed.  It  shows  that  the 
primary  effect  is  to  amplitude  modulate  the  usual  point  source  intensity 
distribution  reflected  from  a target  glint.  With  large  targets,  there  is 
also  an  additive  noise  contribution  from  the  receiver  spatial  intensity 
variations  associated  with  reflections  from  the  area  surrounding  the  tar- 
get focal  spot.  This  analysis  is  also  used  to  explain  the  operations  of 
the  ASPECT-I  and  ASPECT-II  speckle  cancelation  schemes. 

Experimental  data  showing  ASPECT  cancelation  of  simulated  speckle  modu- 
lations is  presented.  The  experiments  show  that  the  laboratory  ASPECT 
system  can  cancel  single-frequency  and  broadband  modulations  within  the 
COAT  dither  band.  ^ 

Previous  work  has  also  shown  that  the  severity  of  COAT  performance 
degradation  due  to  speckle  interference  is  related  to  the  COAT  servo  band- 
width. Generally  speaking,  the  larger  the  bandwidth,  the  more  severe 
the  degradation.  However,  if  the  bandwidth  is  reduced  to  decrease  speckle 
interference,  the  ability  of  the  servo  to  compensate  other  disturbances 
such  as  atmospheric  turbulance  will  also  be  reduced.  Thus,  we 
have  studied  a servo  technique  called  ABAC  that  seeks  to  optimize  the 
COAT  servo  bandwidth  under  arbitrary  conditions  of  speckle  interference 
and  atmospheric  turbulence.  Computer  simulation  results  demonstrate 
ABAC'S  ability  to  achieve  this  goal. 

Analytical  work  which  predicts  the  amount  of  COAT  target  btrehl  ratio 
degradation  for  a specified  but  otherwise  arbitrary  speckle  modulation 
function  is  also  presented.  The  analytical  predictions  are  compared  to 
computer  simulations  and  good  agreement  is  found. 
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SECTION  1 


INTRODUCTION  AND  SUMMARY 


Several  Investigators  have  shown  that  rotating  targets  generate 

speckle  effects,  which  have  the  potciiLi.il  t a seriously  degrade  i , < i 

formance  of  all  adaptive  optics  systems  operating  witli  single  line  laser 
12  3 

illuminators.  ’ ’ Accordingly,  the  primary  objective  ut  this  program 
was  to  investigate  techniques  that  will  give  improved  performance  in  one 
class  of  Lhese  systems  — a multidither  COAT  system  in  ihu  presence  ot 
such  speckle  jamming  i uteri erence . Thus,  we  investigated  several  sui  n 
concepts  and  were  able  to  demonstrate  one  very  promising  concept  that 
can  virtually  eliminate  the  effective  speckle  interference  and  the  tesul- 
tant  COAT  perlormance  degradation.  Inis  concept,  which  we  eall  automa- 
tic speckle  cancellation  techniques  (ASPECT),  refers  Lo  a class  l 
techniques  by  which  speckle  noise  in  a conventional  COAT  receiver  is 
first  estimated  and  then  cancelled.  The  estimation  is  achieved  by  inn 
ducing  a high  frequency  sinusoidal  modulation  the  probe  f requeue i on 
either  a portion  of  the  outgoing  wavefront  (ASPECT  1)  or  on  the  cnLiii 
wavefront  (ASPECT  LI).  The  modulation  envelope  of  this  carrier,  as 
extracted  from  t lie  usual  multidither  detector  output,  provides  tin 
desired  speckle  estimator.  This  estimator  is  used  for  cancellation  in 
one  of  two  ways:  (1.1  the  division  scheme  or  '2)  the  subtraction  schemt 

Basic  systems  configurations  for  ASPECT  1 and  II  and  their  relative 
merits  are  discussed  in  Section  2.  Since  we  concluded  early  in  t lie  pro- 
gram that  ASPECT  11  was  the  preferred  appioa  ii,  it  was  chosen  for  labora 
tory  instrumentation.  ii  these  systems  arc  to  be  useful  in  high  powei 
applications,  ii  is  important  Licit  reasonably  low-cost,  high-power  c-m 
figurations  be  available.  Accordingly,  two  ASPECT  II  conf igurat ions 
suitable  foi  high- powei  applications  arc  alsi>  discussed  in  Secti  n 2. 

To  iindci  -.1  and  how  these  systems  work  requires  a detailed  analysis 
of  tile  speckle  Interference  process  ill  multidither  COAT  systems.  Such  an 
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analysis  is  developed  in  Section  3 and  applied  to  ASPECT  systems  in 
Section  4.  This  analysis  shows  that  the  primary  effect  of  speckle  on 
the  COAT  receiver  output  signal  is  indeed  a multiplicative  noise  that 
can  be  extracted  by  an  ASPECT  I or  II  estimation  system. 

We  recognized  early  in  the  program  that,  with  high  power  systems, 
the  probe  carrier  frequency  would  most  conveniently  be  introduced  via 
phase  modulations  produced  by  moving  mirrors.  Clearly,  these  have 
definite  frequency  limitations  that  might  limit  system  performance. 

More  specifically,  the  lower  the  probe  carrier  freuency  relative  to  the 
(highest)  normal  dither  frequency,  the  higher  the  likelihood  that  the 
probe  modulation  envelope  will  be  contaminated  with  undesired  high- 
frequency  dither  signals  and/or  their  harmonies  as  well  as  the  multipli- 
cative speckle  noise  associated  with  each.  The  ASPECT  frequency  ratio* 
required  to  avoid  such  contamination  is  difficult  to  ascertain  theore- 
tically. This  problem  was  approached  in  two  ways:  (1)  computer 
simulation  and  (2)  laboratory  system  simulation  (as  discussed  in  Sec- 
tions 4 and  5). 

We  had  originally  hoped  that  the  ASPECT  probe  frequency  could  be 
the  second  harmonic  of  the  highest  dither  frequency  (an  ASPECT  frequency 
ratio  of  two)  since  this  would  have  provided  a very  simple,  low-cost 
implementation.  Unfortunately,  the  computer  simulations  indicated  that 
ASPECT  frequency  ratios  of  approximately  six  were  required  for  reason- 
able performance  and  that  ratios  of  ten  or  more  are  required  to  ensure 
the  best  performance.  The  experimental  laboratory  data  confirms  these 
same  conclusions.  In  particular,  although  we  were  able  to  demonstrate 
substantial  improvements  with  the  ASPECT  correction  system  operating, 
we  did  not  see  a total  performance  recovery  as  large  as  that  achieved  in 
the  computer  simulations.  This  was  because  our  maximum  achievable  probe 
carrier  was  100  kHz,  which  was  about  7 times  the  maximum  COAT  dither 


* 

Probe  carrier  divided  by  the  highest  COAT  dither  frequency. 
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frequency  of  15  kHz  (we  used  only  7 low-frequency  dither  channels  in 
the  existing  system  to  improve  the  ASPECT  frequency  ratio),  whereas  10 
times  is  required  for  best  performance. 

We  have  also  conducted  a computer  simulation  study  of  an  automatic 
bandwidth  adaptive  control  (ABAC)  system.  Briefly,  its  operating  princi- 
ple is  as  follows.  The  amount  of  energy  in  the  speckle  interference 
signal  depends  on  the  COAT  servo  bandwidth.  If  we  reduce  the  bandwidth, 
we  reduce  the  interference.  On  the  other  hand,  after  reducing  the  band- 
width, the  COAT  system  may  not  be  able  to  respond  fast  enough  to  compen- 
sate dynamic  phase  distortions  (such  as  atmospheric  turbulence).  The 
ABAC  servo  automatically  seeks  the  optimum  compromise  bandwidth,  the 
bandwidth  that  allows  the  highest  possible  target  irradiance  in  the 
presence  of  these  constraints.  We  have  demonstrated,  using  computer 
simulation,  that  this  concept  does  perform  this  task.  To  implement  ABAC, 
we  dither  the  bandwidth  very  slowly  (about  two  orders  of  magnitude  slower 
than  the  servo  bandwidth)  and  then  synchronously  detect  at  the  dither 
frequency  to  produce  a bandwidth  correction  signal.  Since  the  amount  of 
bandwidth-related  degradation  is  highly  dependent  on  atmospheric  condi- 
tions and  target  scenarios,  the  amount  of  improvement  will  vary 
accordingly. 

In  Section  7,  we  show  that  COAT  target  irradiance  degradation  due 
to  speckle  interference  can  be  predicted  analytically  if  the  power  spec- 
trum of  the  speckle  modulation  function  and  the  COAT  servo  bandwidth 
characteristics  are  known.  The  analytical  predictions  are  compared  to 
computer  simulation  data,  with  good  correlation  between  the  separate 
results. 


SECTION  2 


ASPECT  SYSTEMS 


A.  OPERATING  PRINCIPLES 

ASPECT  is  a series  of  techniques  for  estimating  the  real-time 

jamming  or  interfering  time  modulations  generated  by  a moving  speckle 
1-3 

pattern  and  for  removing  their  effects  (in  real  time)  from  the  adaptive 
optics  error  signals  of  a multidither  system.  This  section  explores 
these  techniques  in  more  detail,  heur istically  describes  their  operation, 
and  considers  their  application  at  high  power  levels.  These  systems 
have  the  following  characteristics  in  common:  (1)  they  share  a single 
receiver  aperature  and  detector  with  the  basic  multidither  system; 

(2)  they  inject  a probing  modulation  or  wavelength  to  extract  the  speckle 
waveform  estimator  via  an  auxiliary  probe  beam  that  shares  the  main 
aperture;  (3)  they  extract  the  modulation  envelope  of  the  probe  frequency 
in  the  detector  output  aS  a real-time  estimator  of  the  speckle  jamming 
signal:  and  (4)  they  divide  this  signal  into  or  subtract  it  from  the 
basic  dither  encoded  signal  out  of  the  detector. 

If  only  a single  element  of  a zonal  multidither  system  is  modulated 
with  the  probe  modulation,  we  call  the  resulting  system  an  ASPECT  I 
system,  while  if  the  entire  (in  the  spatial  sense)  outgoing  beam  is 
modulated,  we  call  the  resulting  system  an  ASPECT  II  system.  Further, 
as  will  be  discussed  in  Section  3.E.  ASPECT  I and  ASPECT  II  estimate 
slightly  different  speckle  noise  components. 

1.  Aspect  I Systems 

The  simplest  approach  is  to  extract  the  ASPECT  probe  frequency  as 
the  second  harmonic  of  the  highest  dither  frequency  in  the  system.  This 
has  the  virtue  that  no  optical  modification  whatever  is  required  of  the 
basic  multidither  outgoing-wave  system  (illustrated  in  Figure  1).  As 
discussed  in  Section  3.E,  it  is  advantageous  to  employ  a central  element 
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in  the  array  for  aspect  encoding.  However,  electronic  modifications 
are  required.  It  is  advantageous  to  raise  the  probe  dither  frequency 
to  the  highest  possible  limit,  operating  at  or  above  a high-frequency 
resonance  of  the  dither  mirror,  or,  if  optical  power  permits,  to  include 
a high-frequency  electrooptic  dither  system  for  this  particular  element. 
Further,  the  probe  dither  amplitude  is  increased  to  boost  the  magnitude 
of  the  second  harmonic  signal.  Specifically,  a boosted  amplitude  tp 

A 

for  the  ASPECT  channel  of 


p.  > 1.8  rad 
A 


(2.1) 


is  recommended  for  this  appiication.  Under  these  conditions,  the 
second  harmonic,  which  is  the  probe  frequency,  is  of  strong  magnitude 
(with  good  convergence  in  this  channel).  Assuming  the  speckle  pattern 
scans  past  the  receive  aperture  and  generates  a speckle  jamming  signal 
Me>(t),  then  a nearly  equivalent  signal  M (t)  also  modulates  the  probe 
frequency.  Thus,  the  high-intensity  components  in  the  detector  output 
current  are  approximately 


S(i  = M^U)  + CMQA(t)  cos(2wAt  + aA)  + ^ Mg^(t)  A^sin^t  +°>n)» 


where  M^^t)  and  M (t)  are  nearly  equal  multiplicative  noise  components, 
S 1 SA 

at  least  for  closely  spaced  glints;  is  the  aspect  probe  frequency, 

the  variables  a and  a,  are  arbitrary  reference  phases;  C is  a constant; 
n A 

and  A is  the  modulation  index  of  channel  n. 
n 

The  envelope  of  the  received  second  harmonic  probe  signal  M (t), 
is  extracted  by  synchronous  detection  and  forms  our  estimation  signal 
M (t),  as  illustrated  in  Figure  2.  This  signal  is  then  divided  into 

bA 

the  basic  common  dither  signal  system  (as  illustrated)  or  substracted  in 


the  alternate  approach.  Since  Mg2(t)  is  the  speckle  noise  signal  of 
major  importance  and  since  MgA(t)  is  not  precisely  equal  to  Mg2(t), 
we  believe  that  these  differences  may  cause  some  reduction  in  perfor- 
mance. Additional  problems  are  discussed  below.  We  use  matched  low-pass 
filters  F^  and  in  Figure  2.  Basically,  the  speckle-estimator  low- 
pass  filter  F^  has  a cutoff  that  must  fall  somewhat  above  the  highest 
dither  frequency  and  below  the  ASPECT  probe  frequency  (2  <uA).  The  filter 
F^  is  included  only  to  balance  the  delay  times  for  the  true  speckle 
signals  and  the  estimator  signal  paths  going  into  the  divider. 

If  the  probe  frequency  were  sufficiently  high  and  no  other  types  of 
noise  were  present,  then  the  speckle  estimator  MgA(t)  and  the  actual 
speckle  multiplier  would  be  substantially  identical.  Unfortunately, 
with  second  harmonic  operation  and  equally  (or  near  equally)  spaced 
dither  signals,  a large  number  of  second  harmonic  signals  from  the  other 
dithers  can  fall  near  the  ASPECT  probe  frequency,  as  illustrated  in 
Figure  3.  In  particular.  Figure  3 illustrates  a "typical"  spectrum  of 
the  dither  channels  with  all  COAT  servo  loops  open,  assuming  the  system 
is  working  well.  The  reference  channel  at*  7 kHz  dominates  in  intensity 
because  of  its  high  modulation  index,  if  we  assume  that  the  dc  component 
is  eliminated  by  high-pass  filtering.  If  the  servo  loops  are  closed,  as 
illustrated  in  Figure  3(b),  the  fundamental  dither  frequency  components, 
which  represent  the  error  signals,  generally  reduce  in  magnitude,  parti- 
cularly at  the  higher  dither  frequency,  in  spite  of  the  speckle  jamming, 
if  the  ASPECT  system  functions  correctly.  In  particular,  the  reference, 
or  probe,  channel  has  a very  large  expected  reduction  at  its  fundamental 
frequency  (good  convergence),  both  because  of  its  high  frequency  and 
because  its  enhanced  modulation  index  makes  it  iess  vulnerable  to  speckle 
interference.  As  the  fundamental  component  reduces  in  amplitude,  the 
second  harmonic  increases  until  it  becomes  the  dominant  spectral  line. 


The  illustrated  frequencies  are  only  examples. 


Figure  3.  Key  voltage  (or  current)  spectra  in  an  ASPECT  I system.  The  indi- 
cated magnitudes  are  discussed  in  the  text.  (a)  Spectrum  of  high- 
pass  filter  output  before  convergence.  (b)  Spectral  output  of 
Figure  22(a)  after  "convergence"  (with  speckle  jamming). 

(c)  Reference  oscillator  spectrum.  (d)  Spectrum  of  reference  syn- 
chronous detector  output,  before  low  pass  filtering. 
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Next  we  examine  the  spectrum  of  the  unfiltered  output  of  the  ASPECT  I 
detector.  For  the  moment,  we  assume  that  our  reference  channel  second- 
harmonic  synchronous  detector  (at  14  kHz)  is  a true  multiplier  and 
produces  an  output  with  the  spectrum  indicated  in  Figure  3(d).  This 
spectrum  is  the  convolution  of  the  reference  oscillator  spectrum  of 
Figure  3(c)  with  the  input  spectrum  of  Figure  3(b).  The  dominant 
spectral  components  are  the  dc  component  and  second  harmonics  at  28  kHz, 
which  are  not  shown  in  Figure  3(d). 

So  far,  the  effect  of  the  speckle  on  the  spectra  in  Figure  3 has 
not  been  illustrated.  It  is  easily  obtained  by  convolving  the  speckle 
spectrum  with  each  of  the  illustrated  spectra.  In  particular,  if  all 
the  spectral  lines  in  Figure  3(d),  (except  the  control,  or  zero  frequency, 

line)  were  vanishingly  small,  then  this  convolutional  process  would 
simply  replicate  the  speckle  spectrum  at  the  detector  output  and  there- 
fore the  temporal  output  of  the  2w  synchronous  detector  would  replicate 

A 

MgA(t).  As  illustrated  in  Figure  3(d),  other  spectral  lines  are  present 
and  contaminate  the  estimation  of  the  speckle  spectra  since  they  couple 
erroneous  spectral  information  via  the  convolution  process.  However, 
the  better  the  estimation,  the  better  the  speckle  can  be  removed  from 
the  lower  frequency  channels,  and  hence  the  better  the  convergence  and 
the  smaller  the  spectral  lines  above  7 kHz  become,  etc.  Thus,  the  effect 
of  the  harmonic  jamming  signals  is  a highly  nonlinear  and  very  complex 
process.  This  problem  is  discussed  in  Sections  4 and  5,  which  describe 
our  computer  simulations  and  our  laboratory  systems  investigations. 

Similar  (but  differing  in  detail)  self-jamming  problems  occur  in 
all  the  ASPECT  systems  diesussed  below  if  the  probe  frequency  is  not 

sufficiently  well  separated  from  the  dither  frequencies.  ♦] 

Using  computer  simulation,  we  investigated  this  jamming  problem 
for  ASPECT  I and  have  found,  under  the  assumptions  of  pure  multiplicative 
noise  as  defined  in  Section  3.C  and  sufficiently  high  probe  frequencies, 
that  the  system  could  be  made  to  perform  well.  However,  because  of  the 
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important  advantages  of  ASPECT  II  (discussed  in  Section  2) , we  elected 
to  devote  the  remainder  of  the  report  to  this  class  of  system. 

2.  ASPECT  II  Implementations 

ASPECT  II  systems  have  demonstrated  the  best  performance  in  the  face 
of  speckle  jamming  of  all  the  systems  considered  in  this  report  and  (so 
far  as  is  known)  elsewhere.  These  are  implemented  by  splitting  off  a 
portion  of  the  whole  beam  (in  the  spatial  sense),  modulating  it  or 
frequency  offsetting  it,  and  recombining  it  with  the  main  beam  ahead  of 
the  wavefront  correction  system  (generally  a deformable  mirror).  Several 
modulation  techniques  potentially  are  applicable,  including  (1)  frequency 
offset  of  the  probe  beam,  (2)  amplitude  modulation  of  the  probe  beam,  and 
(3)  phase  modulation  of  the  probe  beam.  Primarily  because  of  its  high 
recombination  efficiency,  we  selected  the  phase-modulation  approach  for 
initial  discussion.  Further,  since  it  can  employ  mechanically  moving 
mirrors  to  generate  the  probe  modulation,  it  appears  certain  that  it 
can  be  implemented  at  high  power  levels,  without  technology  challenges. 

In  addition  to  the  versatility  in  implementation  techniques  that 
this  approach  permits  in  comparison  with  ASPECT  I,  it  has  the  property 
that  is  directly  estimates  M ^(t)  in  Eq.  2 the  speckle  noise  component 
of  most  importance  in  speckle  jamming.  That  is,  for  ASPECT  II,  M (t) 
is  equal  to  M ^(t)  up  to  a multiplicative  constant.  This  is  discussed 
in  more  detail  in  Section  4. A. 

A typical  low-power  ASPECT  II  configuration  using  phase  modulation 
is  illustrated  in  Figure  4.  In  this  particular  implementation,  which 
closely  resembles  our  laboratory  demonstration  system,  2%  of  the  main 
beam  power  is  split  off  from  the  main  beam,  phase  modulated,  and  then 
recombined  with  the  main  beam  in  a Mark-Zender  interferometer.  A local 
null-seeking  dither  ser'o  holds  the  interferometer  path  lengths  in 
constructive  interference  on  the  recombined  beams*  and  in  destructive 

* 

The  in-phase  operation  of  the  test  beam  serves  two  functions.  It  pro- 
vides the  most  efficient  use  of  the  test  beam  from  a total  delivered 
power  viewpoint,  and  it  also  delivers  the  maximum  output  signal  at  the 
second  harmonic. 
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interference  at  the  local  loop  detector.  Separate  dither  and  corrector 
deformable  mirrors  are  used  and  are  not  simply  interchangeable  in  posi- 
tion or  combinable  without  degrading  system  performance.  In  particular, 
if  the  dither  mirror  is  moved  to  the  position  of  the  deformable  mirror 
(outside  the  interferometer),  the  ASPECT  output  signals  will  include 
full  strength  dither  signals  and  the  division  operation  will  cancel  the 
dither  and  the  speckle  signals.  Although  the  result  of  moving  the 
deformable  mirror  to  the  position  of  the  dither  mirror  is  less  disasterous, 
it  nevertheless  compromises  the  system.  Specifically,  the  probe  beam 
then  becomes  a badly  corrected  beam,  in  the  face  of  propagation  path 
aberrations,  and  will  illuminate  a wider  portion  of  the  target  than  the 
main  beam.  Consequently,  each  beam  can  produce  differing  speckle  fields 
at  the  receiver  aperture  with  differing  associated  amplitude  modulations. 
Thus,  the  ASPECT  cancellation  could  be  less  than  complete,  although 
possibly  still  worthwhile.  Although  both  United  and  Hughes  proposed 
separate  dither  and  corrector  mirrors  in  their  HiCLAS  design  studies, 
there  may  be  future  designs  which  prefer  to  combine  these  functions  in 
one  mirror.  If  so,  the  preferred  location  is  that  of  the  dither  mirror 
(within  the  interferometer)  of  Figure  4.  A low-power  corrector  mirror 
can  be  added  to  the  probe  beam  (in  the  lower  path  of  the  interferometer) 
if  desired  to  correct  the  problem  just  discussed. 

The  receiver  aperture  may  be  off-axis  (as  illustrated),  a shared- 
aperture  annulus,  or  whatever  — it  does  not  matter.  Similarly  the 
receiver  may  be  a simple  "photon-bucket"  collector  or  an  imaging  system 
with  a restricted  field  of  view  — an  "IMPACT"  system.  The  ability  to 
estimate  speckle  is  the  same  in  all  cases  although  the  overall  speckle 
suppression  may  be  influenced  by  these  choices. 

As  in  ASPECT  I,  the  speckle  modulation  envelope  is  extracted  from 
the  second  harmonic  of  the  test  dither  by  synchronous  detection  and 
is  divided  into  (or  subtracted  from)  the  main  error  signals  to  cancel 
the  speckle.  This  is  done  with  an  electronic  system  such  as  described 
for  the  ASPECT  I (See  Figure  2).  For  the  reasons  previously  discussed, 
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it  is  important  to  have  the  probe  frequency  well  above  the  highest 
dither  frequency. 

There  are  some  other  considerations  worth  noting.  For  example, 
there  is  an  optimufti  splitter-recombination  system  for  most  efficient  use 
of  the  reinjected  power.  In  particular  the  i nat i i 

intensity  ratio  shouLd  have  the  same  value  oi  me  intensity  ratios  as 
in  tile  main  and  injected  beams,  with  t lie  majority  ji  the  power  heiag 
transmitted  through  the  splitter  (when  taken  one  he am  at  a time).  in 
this  way,  tin:  two  fields  transmitted  to  t lie  1 ,-al  Loop  detect  a are 
equal  in  magnitude,  and  the  local -loop  system  drives  them  to  an  - ut— oi  ■ 
phase  condition.  this  generates  111  esstnl  iai'c  1 eld  (aua  energy)  null 
in  this  path.  in  other  words,  all  ul  tin.  i n j e ted  powei  it  the  test 

beam  is  transmitted  lu  tlie  target  by  ■ itt  a or  tiiis  inter t ere  a < and 

conservation  oi  energy. 

it  is  tie]  pf  ul  to  define  two  k i ud  s o t t i e 1 J i a 1 i . •.->  in  t li  i • ■ [ 

i lie  t i r s t is  tin.  basic  sp  1 i 1 1 e l t i e ! el  l i 1 1 1 u i a i a 1 . i > t i . . * i 

call  n and  1 1 and  foi  the  case  illustrated  in  Figuie  >: 


nj  - (U.U2)1 

n ( 1 1 . 9H ) 

The  second  factor  is  the  iatie  oi  t lie  field;  in  tin  probe-  i<  itn  I i.c 

main  beam  as  they  exit  Lin  second  splittei  (whie.h  we  assume  ai<  is. ■ i . I . 
as  discussed  above).  these  factors  ate  t lie  quaies  ot  the  q I a.  lot..; 
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The  1 it  Iters  on  the  1st 
mil  ioi  pi  i i u i tin  in  1 1 I 
t rausmi ssiou  ot  e is. 

loe  i i detee  tot  . 


Iiealll  and  t L > ■ a lea., 
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(2.4b) 


For  lossless  splitter,  it  follows  that 


ni  + n2  = _ri 


+ h 2 = 1 


(2.5) 


The  efficiency  of  this  system  is  difficult  to  replicate  with  the 
other  approaches.  However,  these  are  arguments  for  operating  the  test  probe 
beam  at  a second  wavelength  or  at  an  offset  frequency.  The  second  wave- 
length must  not  be  far  separated  from  the  main  beam  wavelength  to  ensure 
that  its  speckle  patterns  replicate  those  of  the  main  beam.  More 
specifically,  if  the  power  levels  permit,  one  may  advantageously  substi- 
tute an  acousto-optic  modulator  for  the  moving  mirror  phase  modulation 
system  of  Figure  4.  The  probe  beam  is  then  up  (or  down)  converted  by  a 
convenient  frequency  (say  1 MHz)  and  recombined  with  the  main  beam,  as 
in  the  Figure  4 system.  In  this  case,  the  up-converted  beam  is  reflected 
and  beats  with  the  reflected  main  beam  in  the  receiver  detector  to  produce 
a 1-MHz  output  — the  probe  frequency.  With  this  approach,  control  of 
the  recombination  path  is  not  required,  which  simplifies  the  system. 

The  main  advantage  is  that  the  probe  frequency  will  be  sufficiently  high 
to  greatly  exceed  any  dither  frequency  range  available  from  present  or 
future  deformable  mirrors.  Thus,  probe  frequency  contamination  problems 
are  avoided. 

B.  HIGH-POWER  IMPLEMENTATION  OF  ASPECT  II  SYSTEMS 

This  section  illustrates  possible  approaches  to  high-power  imple- 
mentations of  ASPECT  systems.  The  objective  is  not  so  much  to  identify 
the  "best"  configuration  as  to  illustrate  the  wide  variety  of  possible 
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approaches.  Thus,  we  select  two  quite  different  implementations.  The 
first  approach,  illustrated  in  Figure  5,  is  a high-power  equivalent  of 
the  basic  system  shown  in  Figure  4 with  diffraction  gratings  on  cooled 
metal  mirrors  substituting  for  the  beam  splitters  and  an  acousto-optic 
frequency  shifter  substituting  for  the  ASPECT  (phase  modulator)  dither 
mirror  (for  higher  bandwidth,  not  higher  power).  To  match  the  inter- 
ferometer beam  sizes  to  botn  the  deformable  mirror  and  acousto-optic 
state-of-the-art  sizes,  we  have  included  a beam  reduction  optics  in  the 
low  power  path.  (We  would  include  such  a reduction  even  with  a dither 
mirror  approach  since  it  reduces  the  size  and  boosts  the  resonance 
frequency  of  this  mirror.)  The  interferometer  as  shown  in  not  a balanced 
path  system  and  is  thus  applicable  only  to  a single  line  laser.  It  can, 
of  course,  be  balanced  by  lengthening  the  upper  path. 

it  is  ackward  (and  expensive)  to  incorporate  a pair  of  high-power 
diffraction  gratings  devoted  exclusively  to  the  ASPECT  function  (after 
the  approach  in  Figure  5).  Fortunately,  one  can  often  override  or 
piggyback  the  ASPECT  system  on  top  of  an  existing  grating  or  gratings 
system.  Thus,  no  new  gratings  are  required.  One  such  approach  is  to 
piggyback  the  system  on  top  of  a MAST  grating,  illustrated  in  Figure  6(a). 
The  normal  function  of  the  MAST  grating  is  to  extract  for  alignment  an 
IR  image  of  the  target  and  a beam  sample  of  the  high-power  beam  at  a 
common  or  near-common  sensor  point.  The  MAST  system  works  best  with 
polarized  laser  beams,  as  illustrated  (because  of  the  diffraction 
sensitivity  to  polarization).  For  our  ASPECT  piggyback  system,  it  is 
also  essential  that  the  laser  be  polarized. 

Figure  6(b)  shows  the  basic  MAST  system  with  an  ASPECT  overlay. 
Briefly,  the  MAST  grating  splits  out  1%  of  the  HEL  beam  in  the  reflected 
order,  encodes  it  with  an  ASPECT  modulation,  rotates  its  plane  of 
polarization.  The  beam  then  reflects  off  the  MAST  grating  with  high 
efficiency  and  is  transmitted  to  the  target  as  illustrated.  The  ASPECT 
encoded  beam  and  the  main  unmodulated  beam  are  both  reflected  from  the 
target  and  beat  together  in  a detector  to  given  an  i.f.  output  that  is 
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Figure  6(a).  , , 

The  basic  MAST  grating  system  extracts  samples  of  the  laser 
beam  and  IR  return  from  the  target  for  purposes  of  angular 
alignment. 


DITHER  AND 
CORRECTOR 
MIRROR 


I 

1 

« Figure  6(b). 

An  ASPECT  II  overlay  on  a MAST  grating  system.  The  anti- 
dither mirror  effectively  removes  the  dither  from  the  low 
power  arm  of  the  interferometer.  An  additional  high  fre- 
quency dither  is  applied  to  all  elements  in  parallel  to 
modulate  the  probe  beam. 
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the  ASPECT  probe  signal.  Since  the  modulated  and  unmodulated  beams 
experience  the  same  Doppler  shift,  no  Doppler  tracking  or  offset  is 
required  even  though  coherent  detection  is  employed  (in  a sense). 

The  modulated  and  unmodulated  beams  exit  the  MAST  with  orthogonal 
polarizations.  If  they  were  to  reflect  from  the  target  and  pass  through 
a collector  aperture  optics  with  these  polarization  states  intact,  they 
would  not  beat  in  the  detector.  Fortunately,  many  target  reflecting 
structures  depolarize  the  returns,  including  (1)  recessed  screws, 

(2)  double  bounce  speculars,  and  (3)  very  rough  areas.  Further,  even 
with  single-pass  (curved-surface)  target  reflectors,  the  reflected 
returns  will  in  general  be  depolarized  as  they  bounce  back  and  forth 
on  the  metal  mirrors  of  a typical  beam  director  or  in  the  receiver 
telescope. 

With  this  approach,  it  is  not  feasible  to  add  a dither  mirror  in 
the  "upper  leg"  of  the  interferometer  because  it  is  not  accessible. 
Therefore,  as  illustrated  in  Figure  6,  the  dither  mirror  and  deformable 
mirror  must  be  located  either  before  or  (better)  after  the  MAST  system 
(outside  the  interferometer).  To  avoid  these  problems,  we  added  a low- 
power  "anti-dither"  mirror  as  a reflecting  element  in  the  ASPECT  encod- 
ing path.  The  same  mirror  can  have  the  ASPECT  dither  added  to  all 
actuators  to  achieve  the  ASPECT  encoding.  Of  course,  for  high  frequency 
encoding,  an  acousto-optic  approach  is  also  applicable  to  this  system. 

We  require  a MAST  grating  that  is  extremely  efficient  in  diffrac- 
tion of  one  polarization  and  reasonably  efficient  in  reflection  for  the 
other  polarization.  Fortunately,  this  is  a natural  tendency  for  most 
designs.  Figure  7 illustrates  a theoretical  prediction  for  a grating 
with  idealized  grooves  (designed  for  9.28  pm)  in  which  all  of  the 
coupling  efficiencies  are  100%.  Measurements  on  real  gratings  indicate 
that  at  least  75%  efficiency  is  available  for  the  second  pass  (rotated 
polarization)  reflection  of  the  modulated  beam  off  of  the  MAST  grating 
and  that  98%  (or  better)  efficiency  is  achievable  for  the  main  beam 
polarization. 
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SECTION  3 


SPECKLE  INTERFERENCE  IN  ZONAL  MULTIDITHER 
COAT  SYSTEMS 


A.  PROBLEM  STATEMENT 

The  analysis  presented  in  this  section  describes  how  speckle  effects 
alter  COAT  receiver  illumination  and  the  resultant  photodetector  current. 
It  shows  that  speckle  is  primarily  a multiplicative  phenomenon.  That  is 
to  say,  the  effect  of  speckle  on  the  COAT  receiver  is  to  multiply  the 
normal  COAT  signal  without  speckle  (i.e.,  a signal  resulting  from  an 
ideal  point  source  target  or  gling)  by  a single  temporal  modulation  func- 
tion that  can  be  determined  from  target  surface  characteristics  and 
target  illumination.  There  is  also  additional  additive  noise  from  the 
light  reflected  off  of  the  target  surface  surrounding  the  glint.  This 
additive  noise  should  be  treated  separately  from  speckle  effects  because 
it  fits  more  into  the  category  of  signal-to-noise  problems  and/or  glint 
hopping  problems.  It  is  always  present  and  is  not  appreciably  altered 
by  including  speckle  effects. 

Consider  a typical  COAT  transmitter  aperture  with  an  area  that  has 
been  divided  into  N elements  or  zones  for  the  purpose  of  adjusting  the 


phase  distribution  of  the  outgoing  optical  field. 


4-6 


Referring  to  the 


geometry  shown  in  Figure  8,  let  U (x')  be  the  incident  complex  scalar 

on  th 

field  on  the  target  plane  from  the  n element,  where  the  vector  x'  is 

the  two-dimensional  position  vector  in  the  target  plane.  The  total 

field  U (x',t)  is  therefore  just  the  sum  of  all  the  element  fields: 
o 


N_ 

V5'-  t>  ‘ I "on^’-O 

n=l 


(3.1) 


The  phase  distribution  of  the  n*1*1  element  field  consists  of  a mean  value 

B (x’.t)  and  a sinusoidal  dither  \p  sin(u>  t+a  ),  where  <1/  is  the  dither 
n n n n n 
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amplitude  in  radians  of  wavelength,  u>n  is  the  dither  frequency,  and  a 
is  an  arbitrary  phase.  Thus,  we  can  express  the  total  field  as 
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V (x* ) e 
on 
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(3.2) 


n=l 


where  V (x')  is  a real-valued  function  representing  the  amplitude  of  the 
th  on 

n element  field. 

The  target  can  be  represented  by  a dimensionless  complex  reflectance 
i d)  i \ 

function  p(x')e  defined  on  the  target  plane,  where  <p  (x1  ) and  p(x') 

are  real-valued  functions  with  actual  values  or  more  conveniently,  sta- 
tistical properties  that  can  be  determined  from  the  characteristics  of 
the  target  surface.  The  field  reflected  from  the  target  plane  is  just 
U^Cx'.t)  multiplied  by  the  reflectance  function.  Assuming  that  the 
receiver  plane  is  in  the  far  field  of  the  target  plane,  the  reflected 
field  can  be  propagated  using  the  Frannhofer  integral.  The  field  at  the 
receiver  plane  is  given  by 


UR(x,t)  = 


ikz  l2z  X " X 
e e 

iXz 


J d2x'  UQ(x',t)  p(x'  )e1<t>(^X  ^ e 1 


. k.  -*• 

- T(x  • x ), 


(3.3) 


where  x is  the  two-dimensional  position  coordinate  on  the  receiver  plane, 
z is  the  distance  from  the  target  to  the  receiver,  X is  the  wavelength, 
and  k is  the  wave  number  of  the  propagated  light.  The  integral  in 
Eq.  3.3  is  easier  to  handle  if  we  assume  that  the  reflectance  function 
is  glint-like: 


D(x' ) = X 6 (x’-  x' ) R , 

i * 


(3.4) 
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where  is  a real  number  with  dimensions  of  area  representing  the 
reflectance  of  the  infinitesimal  area  around  the  point  x^.  This  assump- 
tion is  not  a limitation  because  we  can  approximate  any  continuous 
distribution  this  way  as  accurately  as  we  want.  The  receiver  field 
becomes 


UR(x,t)  = 
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_e e 
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where  4>  (x^)  is  referred  to  simply  as  <f>^.  The  receiver  intensity 

I (x,t)  is  given  by  the  absolute  square  of  Eq.  3.5: 
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Eq.  3.6  contains  all  the  information  needed  to  understand  the 
effects  of  speckle  on  zonal  multimode  COAT  systems.  The  information  is 
very  general,  however,  making  most  of  its  obscure.  Therefore,  some 
special  cases  are  discussed  in  the  following  sections. 

Eq.  3.6  can  also  be  expressed  in  a form  more  compatible  with  Eq . 2.2 
by  separating  the  summations  over  n and  m into  three  terms,  one  associated 


with  a specific  channel  (i.e.,  n = m = v) , one  associated  with  the  double 
summation  over  the  remaining  channels  (reduced  array),  and  one  containing 
the  cross  coupling  terms  between  the  two.  Through  a lengthy  derivation, 
the  COAT  photo-detector  current  can  be  shown  to  be  expressible  as  in 
Eq.  2.2.  Rather  than  do  this,  however,  it  is  easier  to  simply  assume  a 
reduced  array  field,  formed  from  all  channels  except  channel  v,  beating 
against  the  single  element  field  associated  with  channel  v.  This  approach 
is  taken  in  Section  3.E. 

B.  CASE  1 - POINT  SOURCE 

In  much  of  the  early  COAT  work,  before  speckle  was  recognized  as  a 
problem,  the  target  was  idealized  as  a point  source.  With  this  assumption, 
Eq.  3.6  reduces  to  the  same  expression  used  in  those  early  analyses.  Since 
there  is  only  one  glint,  the  summations  over  £ and  p reduce  to  a single  term, 
and  there  is  no  need  to  explicitly  show  the  spatial  dependence  of  various 
functions  of  position  in  the  target  plane.  The  receiver  intensity  becomes 
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i [ 0 sin(oj  t+a  ) - ip  sin  (a)  t + a ) ] , (3.7) 
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where  R is  the  single  glint  reflectance  amplitude.  With  the  additional 

assumptions  that  all  the  dither  amplitudes  iJj  , were  equal,  and  all  the 

field  amplitudes  V , were  equal,  Eq.  3.7  has  been  the  starting  point  in 
om 

much  of  the  earlier  COAT  servo  analyses.  Note  that  the  receiver  intensity 
is  only  a function  of  time  and  that  there  is  no  spatial  variation  in  the 
receiver  plane. 

If  we  expand  the  second  exponential  in  the  summation  of  Eq.  3.7  into 
a Four ier-Bessel  series,  the  receiver  intensity  can  be  expressed  as  the  sum 
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of  a mean  value  C t ) plus  spectral  components  at  all  the  dither  frequencies 

and  their  sum  and  difference  frequencies.  This  is  done  in  Eq.  3.8,  except 
that  only  the  dither  components  and  those  at  twice  the  dither  frequencies 
are  shown.  Higher  harmonics  are  neligible  at  the  small  values  of  ^ 
(typically,  about  0.3  to  0.4  rad)  used  in  most  COAT  systems.  Thus,  the 
receiver  plane  intensity  I^(t)  resulting  from  a single  glint  can  be  expressed 
as : 


where 
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c. 


CASE  2 - SMALL  CLUSTER  OF  CLINTS 


We  now  consider  a small  cluster  of  glints  located  within  an  area 

that  is  small  compared  with  the  diffraction-limited  spot  size  ot  the 

total  COAT  transmitter  aperture.  Therefore,  it  is  reasonable  to  assume 

that  the  illumination  field  amplitude  distribution  V (x')  and  average 

on 

phase  distribution  6 (x',t)  from  a single  element  (with  an  area  much 

n 

smaller  than  the  total  transmitter  aperture)  are  raiily  constant  ovee 
this  target  plane  area.  This  assumption  becomes  more  and  more  accurate 
as  the  number  of  elements  increases  and  the  aperture  area  in  the  trans 
mitter  plane  associated  with  each  decreases.  Thus,  we  ignore  the  spatiaJ 
dependence  of  these  functions  and  write  them  as  constants,  allowing  us 
to  rewrite  Eq.  1.6  as  the  product  of  a glint- like  return  multiplied  by 
a spatiaL  modulation  function  My(x)  due  to  speckle: 
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(1.12) 


wh  :•  is  the  normalization  constant  (Az)  dependent  on  distance 

h*.  w.  eii  the  t i ' get  and  receiver,  and  M ix)  is  just  the  factm  in  it 

s 

t .*  brack  is.  This  expression  is  t lie  same  as  Eq  1.8  in  its  time 


dependence,  but  the  spatial  variation  in  the  receiver  plane  is  no  longer 
constant.  Since  Mg(x)  can  be  seen  by  inspection  to  be  real,  it  can  be 
expressed  as 


M (x)  = 1 ^ R R cos  (K  . x +<t>  ) , 

s i P i P HP  HP 


(3.13) 


where 


K = — (x'  - x'  ) , 

l p z p i 


(3.14) 


(3.13) 


The  spatial  frequency  distribution  in  this  speckle  pattern  is  directly 
proportional  to  the  spacing  between  pairs  of  glints,  and  the  amplitude 
at  each  frequency  is  proportional  to  the  product  of  the  reflectance  at 
each  glint.  The  maximum  spatial  frequency  for  this  case  would  be 

K = — ' (largest  glint  separation). 


Goldf ischer J Goodman, George, ^ and  Erdman  and  Gellert^  have  done 

extensive  studies  which  characterize  speckle  patterns  in  much  greater 

depth  than  these  cursory  comments. 

As  the  target  moves,  the  spatial  modulations  also  move  past  the 

receiver,  causing  temporal  modulation  of  the  glint-like  reflected 

intensity.  Since  the  electric  signal  S from  the  photodetector  receiver 

o 

is  proportional  to  the  average  intensity  incident  on  it,  we  can  write 


S = K M (t)  1 (t), 
o s g 


(3.16) 


P 


where  K is  a proportionality  constant,  and  M^(t)  is  the  aperture 

averaged  temporal  modulation  function  resulting  from  the  motion  of  Ms(x) 

past  the  receiver.  If  this  temporal  modulation  contains  a significant 

amount  of  energy  at  or  near  the  COAT  dither  frequencies,  then  the 

1 2 12 

resulting  adaptive  performance  can  be  degraded.  ’ ’ 

It  is  interesting  to  examine  the  effect  of  aperture  averaging 

associated  with  this  type  of  spatial  frequency  distribution.  This  can 

be  done  very  simply  if  we  asume  a shared  aperture  of  diameter  D,,  for 

T X 

the  COAT  transmitter  and  receiver.  A diffraction-limited  spot  size, 
therefore,  has  a null-to-null  diameter  Dg  that  is  approximately  given  by 


(3.17) 


In  the  case  we  have  considered,  the  largest  glint  separation. 

Ax  = I x'  - x'  | , is  small  compared  to  D . Therefore,  the  spatial  wave- 
P _ 2 tt  s 

lengths  (i.e. , A = — ) associated  with  this  speckle  pattern  are  all 
K 


much  longer  than  the  receiver  diameter: 


(3.18) 


Consequently,  there  will  be  very  little  aperture  averaging. 

As  the  glint  spacing  increases  to  a length  comparable  to  a diffraction- 
limited  spot  diameter,  higher  frequencies  or  shorter  spatial  wavelengths 
appear  in  the  speckle  pattern. 

The  greater  the  spacing  between  two  arbitrary  glints  located  at 

x'  and  x',  the  higher  the  probability  that  the  phases  8 (x ' ) and  8 (x') 

St  p n (t  n p 

associated  with  array  element  n are  appreciably  different.  Since  the 
equality  of  these  phases  was  a crucial  assumption  in  the  derivation  of 
Eq.  3.12,  the  validity  of  that  equation  becomes  suspect.  On  the  other 
hand,  as  the  spacing  increases,  so  does  the  amount  of  aperture  averaging. 


For  spacings  on  the  order  of  half  a spot  diameter  and  larger,  the 
attenuation  due  to  aperture  averaging  is  virtually  total.  In  addition, 
since  the  total  amplitude  at  a given  spatial  frequency  is  proportional 
to  the  number  of  glint  pairs  having  a separation  corresponding  to  that 
frequency  (as  given  by  Eq.  3.14),  the  amplitudes  at  spatial  frequencies 
associated  with  larger  separations  of  glint  pairs  are  smaller  than  with 
those  associated  with  smaller  separations  because  of  the  statistically 
smaller  number  of  possible  pairs.  Thus,  Eq.  3.12  may  still  be  approxi- 
mately correct  even  when  the  small  cluster  of  glints  is  no  longer 
confined  to  an  area  that  is  small  compared  with  a diffraction-limited 
spot.  This  special  case  is  discussed  in  more  detail  in  Section  3.E. 

D.  CASE  3 - SMALL  CLUSTER  OF  STRONG  GLINTS  PLUS  A DIFFUSE  BACKGROUND 

We  now  proceed  to  a very  general  case  in  a hypothetical  COAT  target 
scenario.  Referring  to  Figure  9,  let  ue  represent  the  target  by  a 
reflectance  function  consisting  of  a small  cluster  of  strong  glints  in 
an  area  that  is  small  compared  with  a diffraction-limited  spot  size; 
the  remaining  area  of  the  target  plane  is  covered  by  a uniform  distribu- 
tion of  weaker  glints  having  random  phases  (i.e.,  a diffuse  reflecting 
surface).  To  calculate  receiver  illumination,  it  is  convenient  to 
separate  the  summations  over  the  target  plane  in  Eq.  3.6  into  the  sum 
of  two  separate  summations:  one  over  the  small  cluster  of  strong  glints 
(which  we  identify  by  a plain  summation  sign  as  before)  and  one  over 
the  remaining  glints  (which  we  identify  by  a prime  over  the  summation 
sing).  As  done  in  the  previous  cases,  we  will  drop  the  spatial  refer- 
ences to  the  amplitude  and  phase  distributions  of  the  target  illumination 
functions  on  the  region  where  the  strong  glints  are  clustered.  The 
receiver  illumination  from  Eq.  3.6  becomes 
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The  first  term  in  Eq.  3.19  is  identical  to  Eq.  3.12.  The  remaining  three 
terms  may  be  regarded  as  additive  noise.  This  noise  may  or  may  not  be 
a problem  depending  on  its  relative  strength  compared  to  the  first,  or 
glint,  term.  For  example,  if  the  background  noise  is  comparable  in  mag- 
nitude to  the  glint  return,  the  COAT  system  could  fail  to  converge, 
with  or  without  the  added  complication  of  speckle  interference.  This 
kind  of  problem  would  generally  be  remedied  by  a restricted  field  of 
view  or  IMPACT^  system.  Such  a remedy  is  completely  compatible  with 
the  speckle  problem  because  it  eliminates  the  additive  noise  and  again 
allows  us  to  characterize  the  speckle  interference  as  simply  a multi- 
plicative noise. 

The  distribution  of  spatial  frequencies  in  the  speckle  terms  of 

Eq.  3.19  is  still  dependent  on  spacing  between  individual  glints 

(i.e.,  x-x).  Thus,  the  arguments  presented  in  Section  3.C,  concerning 
P 1 

aperture  averaging  still  hold.  Namely,  the  speckle  pattern  associated 
with  glints  spaced  on  the  order  of  half  a diffraction-limited  array  spot 
diameter  or  larger  will  be  eliminated  by  aperture  averaging.  The  criti- 
cal case,  which  is  considered  in  Section  3.E,  occurs  when  there  are 
strong  glints  that  are  spaced  slightly  closer  than  half  a diffraction- 
limited  spot  diameter  and  thus  are  close  enough  together  to  avoid 
complete  aperture  averaging. 

E.  CASE  4 - M GLINTS  DISTRIBUTED  OVER  THE  WIDTH  OF  AN  ARRAY  PATTERN 

As  shown  by  Eq.  3.6,  the  general  result  for  an  N-element  array 
reflecting  off  of  M scattering  centers  is  quite  complex.  This  section 
considers  a simpler  case,  a system  that  is  operating  near  convergence 
with  small  dither  amplitudes  i|/q,  that  are  the  same  in  each  channel.  We 
consider  the  effect  of  speckle  jamming  on  the  dc  terms  and  on  the  dither 
signals  and  show  that  the  resultant  speckle  multipliers  may  be  signifi- 
cantly different. 

More  specifically,  we  let  all  N-l  elements  except  for  element  n be 
reasonably  well  phased  and  consider  the  beat  between  the  received  field 
from  this  element  and  the  ensemble  receiver  field  from  the  remaining 
N-l  elements,  which  we  call  the  reduced  array.  Such  a viewpoint  describes 


the  key  behavior  of  multidither  systems,  even  if  the  other  elements 

13 

are  not  well  phased  (for  point  targets).  Since  we  take  N to  be 
reasonably  large  (N  > 20) , the  converged  target  plane  field  from  the  N-l 
elements  is  very  close  to  that  of  the  full  array.  For  example,  with 
square  arrays  of  square  elements,  the  reduced-array  target-plane  field 
would  be  (nearly)  a two-dimensional  sinc-like  function  with  amplitude  and 
width  that  are  being  perturbed  in  time  by  a small  percentage  as  a con- 
sequence of  the  small  dithers.  Explicitly,  we  take  the  field  in  the 
target  plane  U^x’.t),  decomposed  into  the  two  field  components  just 
discussed,  which  gives 


U (x' ,t)  = U (x' ) + U (x')  exp  [i  B + ^ sin(u  t)+  — (x"  • x' ) j , 
o ra  e n o vn'z  n 
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where  U (x1)  is  the  reduced  array  target  plane  field, 
ra 
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and  B is  the  average  phase  of  the  reduced  array.*  The  parameter 
_y  t a 

U (x' ) is  the  field  pattern  from  the  n element  and  x"  is  the  coordinate 

6 th  n 

of  the  center  of  the  n element. 

Each  of  the  field  components  in  Eq.  3.20  illuminates  a (possibly) 
large  number  of  scattering  centers,  which  we  have  taken  to  be  "glint- 
like" 6 functions  as  described  in  Eq.  3.4.  Each  component  generates  a 
speckle  field  which  we  call 
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It  is  easy  to  show  that  the  dithers  produce  only  extremely  small 
perturbations  in  the  average  phase. 
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Each  of  these  speckle  patterns  contains  exactly  the  same  spatial  fre- 
quency components  since  the  sources  that  created  each  one  have  exactly 
the  same  locations.  However,  each  frequency  component  is  in  general 
weighted  differently  for  two  reasons:  (1)  the  illumination  fields  from 

the  reduced  array  and  the  n*"*1  element  have  appreciably  different  widths, 
1c  ->•  — v 

and  (2)  the  — (x^  • x')  exponent  in  Eq.  3.20  produces  a nonuniform  phase 
in  the  element  illumination  fields  (except  for  a central  element). 

If  we  sum  these  fields  to  generate  the  total  field  at  the  receiver 
plane,  its  squared  magnitude  is 
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4’e(x) 


(3.24) 


The  dither-modulation  component  becomes  explicit  if  Eq.  3.24  is  rewritten 
in  the  form 
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(3.25) 


The  "dc"  (unmodulated)  and  ac  (modulated)  terms  in  Eq.  3.25  are  given 
to  an  excellent  approximation  as 
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(3.26a) 


Iac  (x>t)  “ 2J0(*o)  |URra(x)l  1 URe(x)  1 sin(*ra  " *e  + Sra"  6n} 

(3.26b) 

Eq.  3.25  corresponds  to  a static  target.  As  the  target  rotates,  fields 
|URra(x)  | , luRe(x)  | , and  4>ra(x)  - <j>e(x)  translate  past  the  receiver 
aperture,  thus  generating  a corresponding  time  varying  output  in  the 
receiver  detector.  Since,  in  general,  each  of  the  coefficients  in 
Eq.  3.25  has  a different  spatial  frequency  distribution,  each  term  pro- 
duces a.  different  temporal  frequency  modulation  in  the  dc  and  ac  detector 
output  terms  of  Eq.  3.26. 

If  the  glint  distribution  is  very  tightly  packed  in  comparison  with 

the  width  of  the  reduced  array  pattern,  then  both  U (x' ) and  U (x' ) 

ra  e 

essentially  will  illuminate  the  glints  uniformly  and  therefore 


lURra«">l  * K 


(3.27) 


In  other  words,  within  a constant,  K,  the  two  multiplicative  speckle 
patterns  are  essentially  equal  and  as  described  in  Section  3.C  there 
exists  a single  speckle  multiplier: 


M (x) 
s 


(3.28) 


On  the  other  hand,  if  some  of  the  major  glints  are  located  on  the 
intermediate  skirts  of  the  reduced  array  pattern,  the  multipliers  may 
be  appreciably  different.  To  explore  these  effects  quantitatively,  we 
will  consider  a specific  example.  The  target  will  consist  of  three 
coplanar  glints,  one  at  the  origin  and  two  symmetrically  displaced  from 
it  by  a distance  x^  in  the  target  plane.  The  glint  magnitudes  are  1, 

P,  and  p,  as  illustrated  in  Figure  10(a) . The  reduced-array  and  element 
field  components  in  the  receiverplane  are  thus 
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(3.29a) 


U ! = P[U  (0)  + pU  (x’>  cos  (— M_+n)], 

Re  i e eg  z 


(3.29b) 


where  P is  a propagation  factor,  and  fi  is  a phase  factor  stemming  from 
the  possible  wavefront  tilt  in  the  element  field  illuminating  the  glints 
(after  Eq.  3.20).  Explicitly, 


k x " 
on 


(3.29c) 


where  x " is  the  central  position  of  the  n element  in  the  array.  From 
on 

here  on,  we  consider  the  most  favorable  case  where  x " is  zero  (a 

on 

central  element).  Thus,  the  spatial  frequency  distribution  in  both 
components  is  a constant  plus  a single  (cosine)  frequency  component. 

The  spatial  (fundamental)  modulation  coefficient  in  these  two 
components  are  def ined  as 


pUraU2 

Ura(0)  ' 


(3.30a) 
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(3.30b) 


The  spatial  (fundamental)  modulation  coefficient  in  the  dc  and  a c 
components  of  Eqs.  3.26a  and  3.26b  is  thus  (approximately)  for  moderate 
to  small  p. 
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The  modulation  coefficients  of  Eqs.  3.31a  and  3.31b  are  not  generally 
the  same.  To  give  a specific  illustration,  consider  a square  array  of 
25  square  elements  such  that  the  reduced  array  patterns  and  element 
patterns  are,  for  good  convergence  on  all  other  elements. 


U (x')  = 4.8H 

ra 


sin[ 56(x' ) J 
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= Hsinc [0 (x' ) ] , 


(3.32a) 


(3.32b) 


where 


0 (x ' ) = (xe  ‘ x'} 


If—  \ 
ra 


(3.32c) 


and  where  x"  is  the  width  of  an  element,  and  x'  is  the  reduced  array 
e ra 

beamwidth*  in  the  target  plane.  The  parameter  H is  another  propagation 
constant,  which  we  do  not  require  in  explicit  form. 

For  this  particular  case,  Eq.  3.31  gives 


Rdc 
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sine  n(-£-) 

+ sine  [3(7^)] 
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(3.33a) 


(3.33b) 


These  spatial  modulation  coefficients  are  transJated  into  temporal 
coefficients  by  the  scan  of  the  spatial  pattern  past  the  aperture. 
However,  the  filtering  action  of  the  aperture  reduces  the  associated 
temporal  modulation  coefficients  (as  discussed  in  Section  3.E  and  3.D), 


which  become 


/ 


Defined  as  half  the  null-to-null  spacing. 
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fr.  = 2W,  psinc(n(T5-)] 


(3.34a) 


x x' 

FRad  = Wa  p {sine  [ n (—-)  ]+  sine  [^(-|— ) ] ) , (3.34b) 

ra  ra 


where  W , the  aperture  weighting  or  filtering  factor,  is 

3 


W = sinc(k  _ x ) , 

a z r 


(3.35) 


and  2x^_  is  width  of  the  receive  aperture  (assuming  a square  aperture). 

For  the  special  (but  typical)  case  where  the  receiver  aperature 
matches  the  transmitting  aperture  (2xr  = 5x")  , Eq.  3.35  becomes 


sine  [n  (-£-£)] 

ra 


(3.36) 


If  we  estimate  the  time-varying  speckle  modulation  from  the  ac  compo- 
nents (as  given  by  an  ASPECT  I estimation  system,  for  example)  and  divide 
through  by  this  estimation,  the  residuai  (fundamental)  modulation  in 
the  ASPECT-corrected  dc  term  is  given  by  the  ratio 


1 + 2W  psinc  [ n (——)  1 sin(w  t) 
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x ' x ' 

1 + p {sine  ( n(^-)]+sinc[-|(^-)]}  sin(u)mt) 
ra  ra 


x x 

1 + Wap{sinc[n(-^-)]-sinc[^  (_2L)]}sin(wn)t) , 
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(3.37) 


where  the  second  form  of  Eq.  3.37  is  valid  for  reasonably  small  values 

of  W p.  For  the  matched  receive  aperture  of  Eq.  3.36,  the  residual 
a 

modulation  becomes 


51 


X 


\esid  p 


ra 


x'  x' 

sinc"[n(^— )]  -sinc[ n (~~ ) ] sinc^-J-)] 


• 5 vx 


ra 


sin(oj  t)  . (3.38) 
m 


As  illustrated  in  Figure  10(b),  there  is  no  residual  modulation  for  small 
glint  spacings  since  the  two  components  in  Eq.  3.38  nearly  cancel; 
similarly,  both  components  become  small  for  widely  spaced  glints.  The 
maximum  problem  occurs  at  a glint-spacing  to  beamwdith  ratio  of  about 

0.6. 

Although  the  analysis  up  to  this  point  has  considered  only  array 
states  near  convergence,  we  have  also  examined  a limited  number  of 
special  cases  for  arrays  in  a general  state  of  convergence  and  have 
found  the  same  general  conclusions  to  be  true.  There  are  the  same 
three  basic  types  of  speckle  interference  previously  discussed.  Second, 
the  multiplicative  speckle  noise  is,  in  general,  different  for  the  dc 
and  ac  components  in  the  detector  output.  Third,  if  the  major  target 
glints  or  reflectivity  elements  are  concentrated  within  a region  which 
is  small  compared  to  the  overall  (converged)  beamwidth,  the  effect  of 
the  differences  in  speckle  noise  is  negligible. 


RESIDUAL  AMPLITUDE  MODULATION 


SECTION  4 


AN  ANALYSIS  OF  ASPECT  II 


The  power  split  in  one  of  the  beam  splitters  is  thus  in  the  ratio  of  ri 2 
to  q , after  the  discussion  in  Section  2. 

* t *■ 

By  "in  phase"  we  mean  the  average  path  length  of  the  two  beams,  excluding 
dither  variations,  differ  by  an  integral  number  of  path  lengths. 


A.  GENERAL  DESCRIPTION 

This  section  describes  how  an  ASPECT  II  system  estimates  the  speckle 
modulation  function  Ms(t)  and  cancels  out  the  speckle  interference.  For 
the  remainder  of  this  section,  we  drop  the  designation  "II"  and  simply 
refer  to  ASPECT.  Consider  the  schematic  of  a COAT/ASPECT  system  in  Fig- 
ure 11.  A small  fraction*  r,  of  the  total  field  from  a COAT  laser  beam 
is  split  off  into  a secondary  beam  and  is  combined  with  the  remaining 
fraction  n^,  in  the  principal  beam,  with  the  combined  beam  proceeding 
through  the  usual  COAT  optics.  Each  element  or  zone  in  the  principal 
beam  is  tagged  with  a dither  frequency,  which  is  then  synchronously 
detected  and  processed  electronically  to  generate  phase  corrections  at 
each  zone  for  the  entire  beam.  The  secondary  beam  is  piston  dithered 

at  a frequency  u>  , which  is  substantially  higher  than  the  highest  COAT 

A 

dither  frequency,  and  with  an  amplitude  of  radians  of  wavelength.  The 
beams  are  then  combined  in  phase  **  and  propagated  to  the  target.  Follow- 
ing the  development  of  the  previous  section,  we  write  the  incident  field 
on  the  target  plane  as  a sum  of  the  N-element  fields  that  make  up  each 


of  these  beams: 
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Figure  11.  Schematic  of  a COAT/ASPECT  optical  system 


The  amplitude  distribution  V (x')  associated  with  each  element  field  is 

on 

the  same  as  in  Section  3,  but  the  phase  distributions,  r (x',t)  and 

n 

f (x' ,t),  have  different  dither  components.  In  the  principal  beam,  repre- 
sented by  the  first  summation  in  Eq.  4.1,  the  phase  distribution  is  the 
same  as  the  typical  zonal  multidither  COAT  system  of  Section  3: 


T = g (x,t)  + ip  sin  (w  t+a  ), 
n n o n n 


(4.2) 


where,  for  simplicity,  the  dither  amplitude  i p is  the  same  for  all  COAT 
channels.  The  phase  distributions  in  the  secondary  beam  have  the  same 
mean  values  as  the  primary  beam,  but  only  one  piston  dither  amplitude  and 
phase,  which  is  the  same  for  each  element.  Thus, 


rn  -I  en(x,t)  + *A  sin  (mAt  + aA), 


(4.3) 


where  is  an  arbitrary  reference  phase.  As  discussed  in  Section  2, 
Eq.  2.5,  we  require  that 


nl  + n2 


1. 


(4.4.) 


Again,  the  reflected  field  on  the  target  plane  is  obtained  by  multiplying 
the  incident  field  by  the  glint-like  reflectance  function  of  Eq.  3.4. 
Propagating  this  field  back  to  the  COAT  receiver  plane  using  the 
Frannhofer  integral  gives  the  field  at  the  receiver  plane: 
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(4.5) 
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The  intensity  distribution  on  the  receiver  plane,  IR(x,t),  is: 
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This  is  similar  to  Eq.  3.9;  it  is  a very  general  expression  that  contains 
all  the  information  needed  to  understand  the  fundamental  operating  princi- 
ples of  a COAT/ASPECT  system.  We  again  consider  the  special  case  of  a 
small  number  of  target  glints  clustered  within  an  area  that  is  small  in 
comparison  with  a diffraction-limited  spot  size  of  the  total  COAT  trans- 
mitter aperture.  Therefore,  following  the  same  procedure  as  we  did  for 
Case  2 in  Section  3,  we  ignore  the  small  spatial  variation  of  the  ampli- 
tude and  phase  distributions  of  the  N-element  fields  over  this  area. 

This  allows  the  receiving  plane  intensity  distribution  to  again  be 
expressed  as  a product  of  a point-source-type  return  multiplied  by  a 
speckle-caused  amplitude  modulation.  Eq.  4.6  becomes 
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Expanding  the  dither  phases  in  Eq.  4.3  into  a Fourier-Bessel  series  allows 
the  entire  dither  spectrum  to  be  expressed  explicitly.  Thus,  Eq.  4.7 
can  be  rewritten  conveniently  as 
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where  M (x)  is  the  same  speckle  modulation  function  as  that  of  Eq.  3.2, 
s 
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The  higher  harmonics  of  Eq.  4.8  have  amplitudes  proportional  to  the 

corresponding  products  of  higher-order  Bessel  functions  of  \b  and  \p  which 

o A 

make  them  negligible  compared  to  the  above  coefficients.  Equation  4.9  per 
mits  estimating  the  additional  loss  in  the  Strehl  ratio  incurred  as  a 
result  of  adding  the  ASPECT  dither.  For  even  a modest  number  of  elements, 
the  second  term  in  Eq.  4.9  dominates  the  first  near-convergence  and  thus 
controls  the  Strehl  ratio.  Under  this  assumption,  the  multiplicative 
Strehl  factor  which  the  ASPECT  dither  introduces  is 


(4  1J) 
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If  we  choose  ipA  near  the  zero  of  j 
A 

Eq.  4.13  becomes 


4 rad)  and  small  r^,  then 
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(4.14) 


Thus  for  2%  of  the  power  split  off  iruo  m,  nSPECT  encodinfa  systi  <s , we 
have  an  SA  = 0.96.  When  the  target  i ' :.c  spatial  speckle  modulation  func- 
tion M (x)  produces  the  same  temper. 1 , function  M (l)  discussed 

in  Section  3.C.  The  signal  from  t . ,.i  ir.ii  tector  becomes 
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where  I'g(t)  is  the  expression  in  a , m Eq.  4.8.  We  u«ate 

the  modulation  function  multiplied  by  the  icient  B'^  ji  .hro- 

nously  detect  the  signal  Sq  ot  the  i 

To  ensure  a high-quality  es  ' m.  pectral  cont  u M ft) 

o 

in  the  region  of  the  COAT  dither  barn;,  ic  lesirable  to  allow  an  ample 
frequency  interval  around  the  earn m • ly , the  carrie  ;h.u  : 1 be 

close  to  an  order  of  magnitude  hi;;  ■ r hr1  . the  highest  COAT  dither 
frequency.  To  have  a strong  signal,  at  h.  :rrier  frequency  it  n also 
desirable  to  have  a rather  large  A n' : • "T  fiber  amplitude,  <|i  As  much 
as  it  rad  of  wavelength  might  be  used  iau  since  the  amplitude  of  ttie  COAT 
dither  signals  are  reduced  by  such  . as  is  the  mean  convt . gence 

level,  IM,  some  optimum  compromise  vaiu  l ultimately  be  selected. 

This  is  readily  apparent  from  Eqs.  jgh  4.14,  which  dttnu  the 

amplitudes  of  the  mean  level  and  ions  at  the  ditliei  t>  quencies 

in  terms  of  the  COAT/ASPECT  sysi 

The  speckle  estimator  sign,.  Id  also  contains  tin  I isic 

dither  signals  (as  a multiplicat i and  thus  the  process  ot  divid- 
ing (or  subtracting)  this  estimu  basic  dither  chaum  1 signals 

has  the  potential  to  wholly  or  pai  : : l the  basic  dithei  inf  notation 
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and  the  speckle  noise  introduced  by  Mg (x) . This  potential  is  quantified 
below.  Combining  the  results  of  Eqs.  4.9  and  4.11  yields  the  amplitude 


modulation  index  for  the  basic  n dither  signal,  which  is  for  large  N, 
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while  the  corresponding  modulation  index  in  the  ASPECT  estimator  output  is 
obtained  from  Eq.  4.10  as 
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If  we  consider  the  case  (which  we  expect  to  be  typical  in  application)  in 
which 


» n2  , 


then  Eq.  4.16a  becomes,  to  an  excellent  approximation. 
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Comparing  Eq  4.16b  to  Eq.  4.16c  shows  that  the  modulation  index  of  the 
estimator  channel  is  nearly  0.5,  the  basic  channel  modulation  index,  for 
all  possible  phase  error  states  of  the  system. 

Consequently,  a system  that  divides  or  subtracts  the  estimated  signal 
will  cancel  about  one-half  of  the  dither  signal  component.  This 
becomes  important  if  thermal  or  shot  noise  components  dominate  spec! 
noise. 

Note,  however,  that  if  one  attempts  to  use  a probe  beam  with  con 

parable  power  to  the  basic  dither  encoded  beam  (that  is,  with  l/< 

then  nearly  complete  dither  cancellation  will  result  and  the  cancellai i . 

of  speckle  noise  will  prove  to  be  a hollow  victory.  For  example,  tak.-  the 

case  of  n,  = n„  and  an  ASPECT  index*  ibA  =2.4  such  that 
12  A 


W - 0. 


Eq.  4.16a  then  becomes 


-4J  (\p  ) J (<Jj  ) V ^ V sin(B  (t)  - 6 (t)) 
1 ro  o o on  om  n m 


N N 


[1  + J2  Op  )]  V V V V cos(B  (t)  - 6 (t)) 
o o ^ ^ on  om  n m 


n=l  m=l 


*This  is  an  excellent  choice  for  most  system  applications. 


(4 , lbd) 


For  low  values  of  t|>  (such  as  those  typically  employed),  Eqs.  4.16  and 
o 

4.16b  are  nearly  identical,  yielding  a system  with  substantial  cancellation. 

For  the  special  case  discussed  in  Section  3.E,  Eq.  4.6  can  again 
be  separated  into  terms  associated  with  a given  element  or  channel  and 
a reduced  array  composed  of  the  remaining  elements  or  channels.  Such  an 
analysis  would  again  show  a different  speckle  multiplier  for  the  ac  and 
dc  terms.  As  discussed  in  Section  2. A,  however,  the  ASPECT-II  scheme 
estimates  the  dc  multiplier,  which  is  the  principal  source  of  degradation. 
Partial  cancellation  of  the  ac  multiplier  will  also  occur  since  the  two 
are  similar  although  not  identical.  In  any  event,  the  ac  multiplier  seems 
to  produce  negligible  performance  degradation,  as  demonstrated  by  computer 
simulation  results.  We  discuss  this  point  in  greater  depth  in  Section  4.B. 
To  show  that  we  indeed  do  estimate  the  dc  multiplier  by  synchronously 
detecting  at  w , we  will  examine  Eq.  4.6  for  the  general  case  of  an 

A. 

arbitrary  distribution  of  glints  in  the  target  plane,  but  this  time  we 
will  ignore  the  COAT  dithers.  This  simplifies  the  result  without  obscuring 
the  central  issue  of  the  speckle  multipliers.  Under  these  conditions, 

Eq.  4.6  becomes 


V**0  = R*Rp  C0S[V*p  + 7 *-(*;-*p)] 

^ A 2 / ~ p 


N N 


(n?  + n~)  ^ 'S  V (x!)V  (x')cos[B  (x't)-B  (x't); 
1 2 2-ion^omp  n *■  mp 


n=l  m=l 


N N 


+ 2^  cosI<PA  sin(mAt  + <*A)  ] £ ^ V^') 


n=l  m=l 


vom(xp)  cos[en(x;,t)  -em(x:,t) 


m p 


(4.17) 


Expanding  cos  [i^A  sin  (wAt  + <*A)  ] into  a Bessel  series  and  retaining  the 
two  lowest  order  terms  gives 
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generate  a detector  current  of  the  form 

So  - Ms2(t)  + CMs2(t)  cos  (2ooAt  + aA).  (4.19) 

Had  we  included  the  dither  terms,  the  general  result  would  be  essentially 
the  same  as  Eq.  2.2  of  Section  2.A.I. 

B.  ESTIMATION  AND  CALCULATION  OF  THE  SPECKLE  MODULATION  FUNCTION 

Consider  the  diagrams  of  the  two  ASPECT  speckle  cancellation  s<  hemes 
in  Figure  12.  The  electric  signal  Sq  from  the  photomultiplier  receiver 
is  conducted  through  two  separate  channels,  one  containing  a ’ow-pass 
filter  and  the  other  a high-pass  filter.  The  low-pass  filter  passes 
all  signals  with  frequencies  in  the  COAT  dither  band  and  below,  atten- 
uating the  ASPECT  carrier  signal,  which  is  significantly  higher  than 
the  highest  COAT  dither.  Therefore,  the  filter  output  S^  can  be  written 
approximately  as 

Sx  = K Mg(t)  (1^  + 1^)  . <4  20) 

This  is  essentially  equivalent  to  the  usual  COAT  signal  (l.e  , without 
ASPECT),  where  1^  contains  the  modulations  at  the  dither  frequencies 
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Figure  12.  Schematic  diagram  of  ASPECT  speckle  cancellation  systems. 
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N 
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n=l 


sin  (oi  t + a ) , 
n n 


(4.21) 


and  I'  is  defined  in  Eq.  4.9. 
n 

The  portion  of  the  signal  emerging  from  the  high-pass  filter,  on 
the  other  hand,  only  contains  signals  at  frequencies  well  beyond  the 
dither  band,  in  particular,  the  signal  at  the  ASPECT  carrier  frequency 
2u)^.  This  signal  is  synchronously  detected  and  low-pass  filtered  by  a 
filter  identical  to  the  one  in  the  other  channel,  giving  the  output  signal 
S^.  This  signal  can  be  expressed  approximately  as 


S2  “ 2 K "s(t)  “s' 


(4.22) 


where  the  factor  of  one-half  comes  from  the  synchronous  detection  process. 

Inspection  of  Eqs.  4.9  and  4.10  shows,  for  a large  number  of  channels, 

that  the  coefficient  B'  is  related  to  the  mean  convergence  level  I'  and 

S M 

the  dither  modulations  1^,  as  follows: 


BS  " 


Ki  ln  + k2  xd’ 


(4.23) 


where  and  are  easily  calculated  constants.  Electronically  dividing 
the  signal  by  in  accordance  with  the  division  scheme  in  Figure  12(a) 
yields 


i.,1  L/i-V|  i 

S2  K1  [ \ h!  ■ 


(4.24) 


where  we  assume  1^  >>  1^.  The  speckle  modulation  function  M^(t)  has  been 
cancelled  by  the  division  process.  By  processing  this  quotient  signal  by 
an  automatic  gain  control  (AGC)  which  normalizes  the  signal  by  its  low- 
frequency  component  (the  first  term  in  Eq.  4.24),  we  get 
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(4.25) 


1 AGC 


1+7  (W 


1 

IM 


which  is  very  similar  to  the  usual  COAT  correction  signal  containing  the 
dither  modulations  normalized  to  the  mean  or  low-frequency  level.  However, 
the  dither  signals  are  reduced  or  illuminated,  depending  on  the  relative 
values  of  and  K^,  as  discussed  in  Section  4. A.  For  typical  cases, 
where  >>  n2>  this  signal  is  reduced  by  a factor  of  0.5.  This  signal 
loss  can  then  be  simply  restored  by  increasing  the  loop  gain  electronically 
An  alternate  way  to  cancel  out  the  speckle  interference  is  shown  in 
the  substration  scheme  of  Figure  12(b).  There,  signals  and  are 
channeled  through  identical  AGCs.  To  give  a simple  illustration  of  how 
this  process  works,  we  separate  the  modulation  function  Mg(t)  into  the 
sum  of  two  components,  one  containing  the  low  frequency  variations  (M^) 
and  the  other  containing  the  remaining  high  frequency  variations  (Mp). 

Then,  using  Eq.  4.20,  we  can  write 

S1  ■ K(MH  + V (Ii  +li>  • K (Vi  + Vd  + Vi  + Vi>-  «'26> 


When  this  signal  is  normalized  by  its  low-frequency  component,  which 
consists  primarily  of  M^I^,  we  8et 


S 


1 


(4.27) 


where  — is  identified  as  the  undesirable  speckle  interference  in  the 

COAT  dither  band.  Similarly,  when  the  signal  S2  is  processed  by  an  AGC, 
we  get 


S 


2 


^2 

K1  XM 


(4.28) 
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Notice  that  this  process  only  cancels  the  speckle  multiplier  of  the  dc 

term,  L in  this  case,  and  leaves  it  as  a factor  in  the  ac  term,  except 
M 

for  normalization  terms.  In  previous  sections,  we  have  maintained  that 
the  ac  multiplier  has  negligible  effects  on  COAT  performance.  The 
viability  of  this  subtraction  scheme  is  crucially  dependent  on  the 
validity  of  this  claim.  Thus,  in  demonstrating  its  good  performance 
(see  Section  4.c),  we  also  verify  that  an  ASPECT-II  system  that  can  only 
cancel  the  dc  multiplier  for  the  case  where  the  two  multiplier  are  differ- 
ent is  a reliable  technique  for  counseling  speckle  interference  for 
all  cases. 

C.  COMPUTER  SIMULATION  RESULTS 

A computer  simulation  program  modeling  an  18-element  COAT  servo  and 

the  ASPECT  speckle  cancellation  electronics  has  been  developed.  It  is  a 

modification  of  a code  used  in  previous  studies  of  multidither  COAT 
12  12 

systems.  ’ ’ Figure  13,  a schematic  of  the  modified  program,  includes 
all  of  the  electronic  servo  parameters,  but  no  propagation  effects.  The 
receiver  is  ideal  in  the  sense  that  it  detects  the  return  from  a point 
source  (i.e.,  a target  glint).  Speckle  is  then  introduced  as  an  ampli- 
tude modulation  of  the  normal  receiver  signal.  The  speckle  modulations 
themselves  are  calculated  functions  of  time,  generated  by  computer  at  the 
General  Research  Corporation  and  discussed  in  more  detail  in  Ref.  1.  They 
represent  the  aperture  averaged  intensity  fluctuations  at  a COAT  receiver 
caused  by  the  motion  of  a speckle  pattern  from  a rotating  target  sphere. 
The  surface  of  the  sphere  was  assumed  to  be  optically  rough,  thus  causing 
a speckle  pattern  with  a unit  contrast  ratio,  or,  equivalently,  very 
severe  modulation  depths.  Three  different  target  rotation  rates  were 
considered:  0.4  rad/sec  (which  contained  most  of  its  spectral  energy 

below  the  COAT  dither  band),  2.0  rad/sec  (which  concentrated  most  of  its 
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Figure  13.  Schematic  diagram  of  COAT/ASPECT  computer  simulation  model. 
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Figure  14.  Ideal  COAT/ASPECT  performance  — no  speckle. 


Figure  15.  Speckle  degraded  COAT  performance  ASPECT  open  loop  target  rotation  rate  0.4  rad/sec. 
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Figure  18.  Restored  COAT  performance  ASPECT  closed  loop  target  rotation  rate  0.4  rad/sec 


7772-14 


O 


O CM  rt  <D  00  O 

* CO  oi  «"*’  O 

Nouvinaow  3ixo3dS  lvruov 


80 


Figure  21.  Speckle  modulation  function  used  in  computer  simulations.  Target  rotation  rate,  0.4  rad/sec 
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23.  COAT/ASPECT  performance  with  speckle  ASPECT  dither  frequency,  200  kHz  division 
scheme.  Target  rotation  rate  0.4  rad/sec. 


Figure  24.  COAT/ASPECT  performance  with  speckle  ASPECT  dither  frequency,  200  kHz.  Subtraction 
scheme  target  rotation  rate  0.4  rad/sec. 
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Figure  26.  COAT/ASPECT  performance  with  speckle  ASPECT  dither  frequency,  100  kHz.  Subtraction 
scheme  target  rotation  rate  0.4  rad/sec. 


of  both  schemes  at  100  kHz.  The  average  Strehl  ratio  for  the  division 
scheme  is  noticeably  lower,  but  for  the  subtraction  scheme  it  remains 
quite  close  to  ideal.  Simulations  using  an  ASPECT  dither  of  50  kHz  were 
also  run.  They  showed  a further  degradation  of  performance  of  both 
schemes.  In  fact,  at  this  frequency  the  ASPECT  system  actually  caused 
more  degradation  than  was  originally  present  from  the  speckle  modulations 
alone.  The  explanation  for  this  behavior  is  presented  below. 

Consider  the  simplified  hypothetical  power  spectrum  (based  on 
Eq.  4.8  and  shown  in  Figure  27(a))  of  a COAT  receiver  signal  coming  from 
a glint-like  target  return.  There  is  a strong  mean  level  at  zero  fre- 
quency, discrete  dither  signals,  and  additional  signals  at  twice  the 
dither  frequencies  displayed  on  an  arbitrary  frequency  scale.  There 
must,  of  course,  be  signals  at  all  possible  sums  and  differences  of  dither 
frequencies,  but  their  amplitudes  are  very  small  compared  to  those  shown. 
In  addition  to  this,  there  is  a discrete  signal  at  twice  the  ASPECT 
dither  frequency.  Figure  27(b)  shows  a hypothetical  speckle  modulation 
spectrum.  Since  the  receiver  signal  is  the  product  of  the  speckle  modu- 
lation function  and  the  usual  glint-like  return,  its  power  spectrum  must, 
by  the  convolution  theorem,  be  just  the  convolution  of  the  former  two 
spectra.  This  is  shown  in  Figure  27(c).  Notice  that  as  long  as  the 
value  of  ASPECT  carrier  freuency  2to^  is  significantly  beyond  the  discrete 
signals  associated  with  the  glint  return,  the  convolution  process  will 
reproduce  the  speckle  modulation  spectrum  up  to  a multiplicative  coeffi- 
cient centered  around  the  carrier  frequency.  The  synchronous  detection 
process  then  produces  a signal  that  is  instantaneously  directly  propor- 
tional to  the  speckle  modulation  function  and  also  produces  high  frequency 
noise,  which  can  easily  be  filtered.  As  the  ASPECT  dither  frequency  and 
therefore  the  carrier  frquency  are  reduced,  there  will  be  more  and  more 
overlap  between  the  frequency-shifted  speckle  modulation  spectrum  centered 
at  the  carrier  frequency  and  the  remaining  lower  frequency  convolved 
spectrum.  Therefore,  the  synchronous  detection  process  will  detect  the 
erroneous  overlap  spectrum  as  well  as  the  desired  one,  leading  to  a poorer 
estimate  of  the  speckle  modulation  function. 
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This  also  suggests  why  the  subtraction  scheme  works  better  than  the 
division  scheme.  A careful  comparison  of  the  two  schemes  will  show  that, 
if  a perfect  estimate  of  the  speckle  modulation  function  were  available, 
the  division  scheme  would  cancel  it  completely  while  the  subtraction 
scheme  would  leave  residual  uncancelled  noise.  On  the  other  hand,  when 
the  estimate  is  not  perfect,  there  will  be  some  residual  error  using 
either  method.  Since  the  subtraction  scheme  is  far  less  sensitive  to 
error,  it  performs  better  under  such  circumstances. 

To  allow  for  a direct  comparison  between  the  experimental  results  and 
the  computer  simulation  runs,  we  adjusted  the  input  parameters  to  the  simula- 
tion code  to  represent  the  actual  laboratory  system  as  closely  as  possible. 
The  number  of  COAT  channels  was  reduced  from  18  to  7 having  the  same  dither 
frequencies  as  the  laboratory  model  described  in  Section  5. A.  The  ASPECT 
dither  frequency  was  set  at  50  kHz,  and  the  corner  frequencies  of  the  various 
high-  and  low-pass  filters  in  the  ASPECT  electronics  where  adjusted  to  the 
appropriate  values.  The  division  scheme  of  Figure  12(a)  was  used,  since  that 
was  also  the  experimental  configuration.  The  high-pass  filter  consisted  of 
four  poles  with  corner  frequencies  of  19.3  kHz  each,  while  the  two  matching 
low  pass  filters  had  two  poles  each,  and  each  pole  had  a corner  frequency  of 
15.0  kHz.  The  COAT  dither  amplitude  was  0.349  rad  of  wavelength  in  each 
channel,  while  the  ASPECT  piston  dither  amplitude  was  one-half  a wave.  A 
fraction  02  = 0.13  of  the  optical  field  was  used  in  the  ASPECT  beam;  the 
remainder,  0^  = 0.87,  was  channeled  through  the  usual  COAT  optics.  Simula- 
tions were  run  using  a single  sine  wave  speckle  modulation  function  at  two 
different  frequencies,  one  at  6.5  kHz  and  one  at  13.4  kHz.  In  each  case,  the 
amplitude  of  the  sine  wave  was  0.8  and  the  mean  value  was  1.0.  Both  of  these 
modulations  significantly  degraded  the  COAT  convergence  level,  or  target 
Strehl  ratio.  In  Figure  28  we  show  the  ideal,  no  speckle,  COAT  target  Strehl 
radio  for  this  configuration.  In  Figures  29  and  30  we  show  how  this  ideal 
performance  is  degraded  by  the  respective  sine  wave  modulations.  With  the 
ASPECT  loop  opened,  the  average  Strehl  ratio  at  6.5-kHz  modulation  is 
reduced  to  0.36  of  the  ideal  value;  at  13.4  kHz,  it  is  reduced  to  0.51  of 
the  ideal  value.  These  two  numbers  can  be  directly  compared  to  the  experi- 
mental results  in  Section  5.B,  where  we  observed  a degradation  to  0.39  of 
the  ideal  value  at  both  of  these  frequencies. 
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Figures  31  and  32  show  the  ASPECT  closed-loop  results  for  the  same  two 
modulation  functions.  At  6.5  kHz,  performance  is  almost  fully  restored, 
reaching  0.93  of  the  ideal  level.  At  13.4  kHz,  the  restoration  is  to 
0.88  of  ideal.  This  can  be  compared  to  0.91  of  ideal  at  6.5  kHz,  and 
0.65  of  ideal  at  13.4  kHz,  which  are  the  corresponding  experimental  values 
from  Section  5.B. 

D.  CONCLUSIONS 

The  computer  simulation  results  clearly  demonstrate  that  it  is 
possible  to  restore  speckle-degraded  performance  of  multidither  COAT  sys- 
tems using  an  ASPECT  cancellation  scheme.  High-quality  estimates  of  the 
speckle  modulation  function  can  be  generated  if  the  ASPECT  carrier  fre- 
quency is  substantially  higher  than  the  highest  COAT  dither  frequency. 

If  the  ASPECT  dither  frequency  is  too  close  to  the  COAT  dither  band,  then 
poorer  estimates  and  consequently  poorer  cancellation  will  be  achieved. 

As  a general  rule,  our  data  indicate  that  the  frequency  interval  between 
the  highest  COAT  dither  and  the  ASPECT  carrier  should  be  no  less  than 
6 and  prefereably  about  10  times  the  highest  COAT  dither  frequency.  This 
will  ensure  that  the  available  bandwidth  around  the  carrier  will  be  broad 
enough  to  faithfully  reproduce  the  true  spectrum  of  the  speckle  modula- 
tion function  up  through  the  entire  COAT  dither  band.  Beyond  that  fre- 
quency range  it  is  no  longer  necessary  to  cancel  any  of  the  speckle  modula- 
tions because  they  do  not  produce  any  significant  interference.  The 
spectral  range  of  the  speckle  modulation  function  is  also  important.  But 
since  it  depend  on  such  parameters  as  target  motion,  geometry,  and  sur- 
face characteristics,  it  cannot  be  anticipated  a priori. 

There  are  several  other  important  system  parameters  besides  dither 
frequency  that  will  influence  the  quality  of  the  ASPECT  estimates.  The 
amount  of  energy  in  the  secondary  beam  and  the  amplitude  of  the  ASPECT 
dither  are  both  very  important  in  determining  the  strength  of  the  carrier 
signal.  To  achieve  a strong,  easily  detectable  carrier  it  is  desirable 
to  dither  with  an  amplitude  of  about  half  a wave.  If  the  percentage 
power  in  the  probe  beam  is  large,  it  would  have  the  effect  of  reducing 
the  mean  Strehl  ratio  or  convergence  level  as  well  as  the  strength  of 
the  normal  COAT  error-correction  signals  at  the  dither  frequencies. 
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31.  ASPECT  restored  COAT  performance  with  6.5  kHz  sinewave  modulation 


ASPECT  restored  COAT  performance  with  13.4  kHz  sinewave  modulation. 


Similarly,  increasing  the  energy  in  the  ASPECT  probe  beam  would  increase 
the  strength  of  the  carrier  at  the  expense  of  the  convergence  level  and 
dither  signals.  Since  we  were  primarily  interested  in  demonstrating  the 
ASPECT  concept  and  supporting  the  experimental  effort,  no  attempt  was 
made  to  optimize  these  parameters. 

We  also  found  that  the  subtraction  scheme  for  canceling  speckle 
interference  is  generally  better  than  the  division  scheme.  The  division 
scheme  is  a mathematically  more  precise  scheme  in  the  sense  that,  given 
a perfect  estimate  of  the  speckle  modulation,  it  would  achieve  perfect 
cancellation.  The  subtraction  scheme  on  the  other  hand  would  leave 
residual  uncanceled  interference,  even  with  a perfect  estimate.  But 
since  the  estimates  are  never  perfect,  the  subtraction  scheme  usually 
performs  better  because  it  is  less  sensitive  to  errors. 

The  success  of  the  subtraction  scheme  also  demonstrates  that,  for 
the  special  case  discussed  in  Section  3.E,  in  which  the  speckle  inter- 
ference cannot  be  characterized  by  a single  multiplier,  cancellation  of 
the  dc  multiplier  is  all  that  is  really  ncessary  because  the  ac  multiplit- 
will  produce  negligible  interference  in  the  COAT  convergence  process. 
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SECTION  5 


ASPECT  EXPERIMENTAL  STUDIES 

A.  EXPERIMENTAL  HARDWARE  AND  ELECTRONICS 

The  ASPECT  system  used  in  our  experimental  studies  is  an  add-on 
system  intended  to  "piggyback"  the  present  laboratory  COAT  electronics 
thereby  forming  a complete  adaptive  optics  speckle  cancelling  system.  In 
the  optical  arrangement,  the  existing  optics  and  light  paths  surrounding 
the  operation  of  the  present  COAT  deformable  mirror  remain  unchanged 
except  for  the  inclusion  of  two  beam  splitters,  which  are  used  in  the 
ASPECT  optical  channel  (see  Figure  33  for  the  ASPECT  optical  layout). 

The  light  source  for  the  experiment  was  an  argon  laser  tuned  to 

O 

4880  A and  equipped  with  a rotatable  1/2  A plate  and  etalon.  The  beam 
coming  from  the  laser  goes  to  a spatial  filter  and  lens  (LI)  where  it  is 
expanded  and  collimated  and  directed  to  the  first  beam  splitter  (BS1). 
This  is  the  point  of  origin  for  the  ASPECT  probe  path,  which  is  marked  by 
the  broken  line  in  Figure  33.  The  alternate  COAT  path  is  marked  in 
Figure  33  by  a solid  line.  Starting  at  the  beam  splitter  and  following 
along  that  path,  we  first  come  to  lenses  L2  and  L3.  These  lenses  form 
an  expanding  telescope  to  fill  the  deformable  mirror  with  collimated 
light.  After  leaving  the  telescope,  the  beam  passes  through  another 
telescope,  consisting  of  lenses  L3  and  L4,  which  reduces  the  beam  to  its 
original  size.  After  the  light  passes  through  lens  L4,  it  is  again 
collimated  and  directed  to  the  second  beam  splitter  BS2,  after  which  it 
enters  lenses  L5  and  L6,  forming  the  final  telescope.  After  leaving 
lens  L6,  the  beam  illuminates  the  target  (a  small  piece  of  Scotchlite  in 
a field  of  black  velvet),  and  the  spectral  return  from  Scotchlite  is 
received  by  photomultiplier  (PMT1). 

The  ASPECT  light  path  begins  at  beam  splitter  BS1,  where  the  light 
is  directed  to  mirror  M3  and  into  lens  L7  (the  first  lens  of  a 1 to  1 
telescope).  Since  the  dither  actuator  mirror  (PZT1)  is  smaller  than  the 
beam  entering  lens  L7,  it  is  positioned  just  before  the  focal  point  of 


97 


r*r 


r 


7772-58 


Figure  33.  ASPECT  optical  layout. 
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that  lens  and  directed  with  minimum  angle  toward  tne  corrector  mirror 
(PZT2),  which  is  positioned  just  after  the  local  point  of  lens  L7.  The 
light  is  then  propagated  through  Ld,  the  ouLput  lens  of  the  telescope. 
Actuator  mirror  (PZl'l)  encodes  the  ASPECT  beam  with  a 50-kliz  dither 
signal,  which  is  then  synchronousiy  detected  in  the  same  manner  as  the 
usual  COAT  dither  signals.  The  resultant  signal  from  the  synchronous 
detection  process  then  becomes  the  error  signal  which  is  applied  to  the 
path  balancing  corrector  mirror,  PZT2,  which  performs  the  path  length 
adjustments,  ensuring  that  the  ASPECT  beam  is  properly  phased  relative 
to  the  COAT  beam  at  the  target  glint.  After  leaving  lens  L,8 , the  encoded 
ASPECT  beam  passes  through  an  aperture,  where  its  size  is  adjusted  to 
accommodate  a proper  overlay  with  the  COAT  beam  at  the  beam  splitter 
(BS2).  In  the  far  field,  mirrors  (M4)  and  (Mb)  provide  the  position  and 
angle  adjustments  needed  to  control  the  interference  pattern  between  the 
two  beams.  Before  the  beam  leaving  mirror  (M5)  reaches  beam  splitter 
(BS2),  it  passes  through  attenuator  (ATI),  where  the  intensity  ratio 
between  the  two  beams  may  be  adjusted.  In  the  experiments  described  in 
Section  5,B,  we  typically  measured  power  Levels  of  approximately  Q.1S  mW 
in  the  ASPECT  beam  and  7.0  uiW  in  the  COAT  beam.  These  measurements  were 
made  following  lens  L5  in  the  optical  layout,  one  beam  being  totally 
blocked  while  the  power  in  the  other  was  measured. 

As  shown  in  Figure  34,  the  ASPECT  elect  rani • s system  consists  f 
three  separate  components:  (1)  a local  loop  subsystem,  (2)  a speckle 
simulation  subsystem,  and  (3)  a speckle  detection  loop.  The  local  loop 
subsystem  is  essentially  a copy  of  any  one  of  tb*  COAT  channels.  With 
the  exception  of  it3  higher  dither  frequency,  all  other  parameters  remain 
the  same.  A 50-kHz  master  oscillator  supplies  the  input  to  the  local 
loop's  dither  driver,  which  in  turn  tags  the.  ASPECT  optical  beam  with  a 
± \/2  phase  dither  via  the  ASPECT  dither  mirror.  This  oscillator  also 
supplies  the  input  to  the  adjustable  phase  shifter,  which  compensates 
for  the  phase  delay  in  the  dither  driver  eiei  ■'  tonics,  PZT  actuators,  and 
preamps.  The  utput  of  the  phase  shifter  i •,  the  reference  input  to  the 
local  loop's  synchronous  detector,  and  the  ; ■ i L of  tin  detector  drives 
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a 10-Hz  low-pass  filter  to  provide  a correction  signal  for  the  ASPECT 
corrector  mirror  (PZT2  of  Figure  33).  The  signal  input  the  local  loop 
synchronous  detector  originates  at  the  target  PMT.  Prior  to  synchronous 
detection,  this  signal  is  amplified  and  high-pass  filtered  to  remove  the 
undesired  signals  below  50  kHz.  Two  controls  are  provided  in  the  local 
loop  system:  a preamp  gain  control  and  a correction  signal  gain  control. 

The  speckle  simulator  shown  in  Figure  34  represents  the  multiplica- 
tive speckle  interference  discussed  in  varying  detail  in  Sections  2,  3, 
and  4.  It  was  necessary  to  use  an  artificial  multiplicative  noise  source 
rather  than  a true  speckle  pattern  because  of  a signal  to  noise  problem 
in  the  target  PMT  and  associated  preamps.  Even  at  full  laser  power,  the 
receiver  aperture  of  the  PMT  had  to  be  quite  large  (an  order  of  magnitude 
larger  than  the  transmitter  aperture  mm)  in  order  to  operate  in  the  shot- 
noise-limit  range  of  the  photomultiplier.  Although  this  large  aperture 
provided  a very  good  signal  to  the  COAT/ASPECT  electronics,  it  also  pro- 
duced enough  aperture  averaging  of  the  speckle  pattern  to  preclude  any 
speckle  interference.  Therefore,  we  chose  the  following  scheme  for  elec- 
tronically generating  a multiplicative  noise  that  would  have  a power 
spectrum  similar  to  that  observed  from  the  backscatter  of  a rotating, 
optically  rough  sphere. 

We  used  noise  diode  as  a broad-band  noise  source  and  an  adjustible 
amplifier  to  give  amplitude  control.  This  noise  has  a zero  mean  value  and 
a flat  spectrum  well  beyond  the  highest  COAT  dither  frequency.  To 
generate  a strictly  positive  speckle  modulation  function,  the  noise  was 
then  passed  through  an  analog  multiplier  to  produce  the  square  of  the 
original  noise.  This  was  followed  by  a selectable  low-pass  filter,  which 
allowed  for  modeling  the  frequency  roll  off  characteristics  always  present 
in  any  speckle  pattern,  and  in  particular  those  patterns  associated  with  a 
rotating  target  sphere. 

Since  a true  speckle  noise  always  contains  effects  of  apreture 
averaging,  which  prevent  it  from  ever  reaching  zero,  we  also  introduce  a 
small  dc  offset  in  the  simulated  noise  to  ensure  a consistently  positive 
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signal.  Figure  35(a)  shows  a scope  trace  of  a typical  noise  signal  used 
in  the  laboratory  experiments;  Figure  35(b)  shows  its  spectrum  on  a spec- 
trum analyzer. 

In  addition  to  the  noise  described  above,  the  speckle  simulator  could 
also  inject  a single  sine  wave  with  a dc  offset,  or  just  a dc  signal 
alone.  The  latter  was  used  for  comparative  purposes  to  show  how  the 
systems  operated  without  speckle  interference. 

In  the  speckle  detection  loop,  the  output  of  the  target  PMT  is  first 
multiplied  by  the  noise,  sine  wave,  or  dc  signal.  It  then  takes  two  paths, 
one  to  the  15-kHz*  low-pass  filter  in  the  numerator  leg  of  the  speckle 
divider,  and  the  other  through  a 30-kHz*  high-pass  filter  where  it  becomes 
the  signal  input  of  the  speckle  synchronous  detector.  With  the  inter- 
ferometer loop  converged,  this  signal  is  essentially  a 100-kHz  carrier 
whose  modulation  envelope  is  the  speckle  noise.  The  numerator  path  is 
almost  identical  to  the  usual  COAT  signal  processing  without  ASPECT,  the 
only  difference  being  the  low-pass  filter  that  was  inserted  to  match  or 
duplicate  the  signal  distortion  introduced  by  an  identical  filter  in  the 
denominator  leg.  This  assures  that  the  division  process  will  be  time 
synchronous . 

The  reference  signal  to  the  ASPECT  speckle  synchronous  detector 
originates  at  the  50  kHz  master  oscillator.  It  is  first  passed  through  a 
0 to  360°  phase  shifter,  which  adjusts  for  the  total  phase  lag  in  the  speckle 
detection  loop.  The  output  from  this  phase  shifter  goes  directly  into  a 
multiplier  configured  as  a frequency  doubler,  creating  a 100-kHz  sine  wave 
which  is  the  actual  reference  input.  On  synchronous  detection,  the  multi- 
plicative noise  is  recovered  and  additional  high-frequency  product  signals 
are  introduced.  The  matching  15  kHz  low-pass  filter  removes  these  high- 
frequency  products,  leaving  only  the  multiplicative  speckle  noise,  which 
then  serves  as  the  denominator  input  to  the  speckle  divider.  Depending  on 
the  fidelity  of  the  synchronously  detected  and  filtered  speckle  noise 


1 


* 


3 dB  point. 


* 


102 


0.5  MSEC 


(a)  TYPICAL  NOISE  SIGNAL  AS  A FUNCTION  OF  TIME 


(b)  SPECTRUM  OF  (a).  HORIZONTAL  SCALE  IS  2 kHz/DIV 
VERTICAl  SCALE  LINEARLY  NORMALIZED 


estimate  with  respect  to  the  original  speckle  noise  coming  from  the  speckle 
simulator,  we  can  achieve  partial  or  complete  cancellation  by  the  division 
process.  Since  introducing  the  matching  low-pass  filters  generated  a phase 
lag  in  the  output  signal  from  the  speckle  detection  loop,  the  relative  phase 
in  the  synchronous  detector  oscillators  of  the  higher  frequency  channgles  of 
the  COAT  system  had  to  be  adjusted  accordingly. 

B.  EXPERIMENTAL  RESULTS 

When  the  ASPECT  system  was  first  attached  to  the  existing  COAT  sys- 
tem, and  before  any  attempt  was  made  to  get  experimental  data,  the 
linitations  of  the  hardware  were  made  evident  by  limited  frequency  response 
in  the  available  amplifiers,  drivers,  and  actuators.  Because  of  these 
limitations,  50  kHz  was  the  maximum  achievable  dither  frequency  for  the 
ASPECT  dither  mirror.  Consequently,  the  maximum  available  probe  carrier 
frequency  used  in  the  speckle  detection  loop  was  100  kHz.  According  to 
the  computer  simulation  results  discussed  in  Section  4.C,  this  would  be 
too  close  to  the  highest  COAT  dither  (32  kHz)  to  provide  an  adequate 
estimate  of  the  speckle  modulation  function.  A minimum  frequency  ratio* 
of  about  6 or  7 was  required;  the  actual  value  was  3.13. 

To  verify  these  conclusions  experimentally,  several  attempts  to 
produce  reasonable  estimates  of  single  sine  wave  speckle  modulations  in 
the  COAT  dither  band  were  made.  Even  though  the  signal  from  the  speckle 
simulator  always  had  positive  polarity,  the  estimated  signals  would  often 
drop  to  zero  and/or  change  polarity.  This  type  of  signal  could  obviously 

not  be  used  as  the  denominator  input  to  the  divider  without  introducing  i 

large  errors.  The  estimator  output  was  also  very  noisy,  containing  large 
amounts  of  higher  frequency  modulations  in  addition  to  the  particular 
sine  wave  being  replicated,  in  spite  of  the  low-pass  filtering  following 

. 

the  synchronous  detector  that  generates  the  estimates.  An  attempt  to 


Probe  carrier  frequency  divided  by  the  highest  COAT  dither  frequency. 
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close  the  loop  revealed  unstable  r *' luetuations  in  the  receiver 

PMT  signal. 

The  separation  of  the  100-kHx  pro':  ri  i.  i t requency  from  the  highest 

COAT  dither  (32-kHz)  was  not  lur,  . t <‘n orate  a good  speckle  esti- 

mation signal.  To  achieve  the  m . ess.i  tn  iucn.  separation,  the  number  of 
active  COAT  channels  had  to  be  rei.u.  . . t to  7 with  the  highest  dither 

frequency  being  15  kHz.*  This  gave  in  ratio  of  b.67,  which  was  in 
the  range  that  should  allow  a reasonable  < -tiniate  and  appreciable  cancella- 
tion. However,  the  computer  simulation  data  till  suggested  that  complete 
cancellation  was  not  to  be  expected. 


The  above  modifications  required  lxc.i  tv.  adiustmeuts  in  the  COAT 
servo  parameters  that  are  too  numerous  t<  let  i:  nere.**  Suffice  it  to 
say  that  the  COAT/ASPECT  system,  inch,  it  . rt : il  paths  and  electronics, 

was  once  again  optimized  to  the  ; ’ at  w . ■ t -.vstems  performed  well 
without  speckle  interference  aru  • >.  , ,t  iuatnr  could  be  monitored 

and  observed.  On  introducing  i totally  a 'live  polarity  and  single 
sine  wave  plus  dc  offset  distuihance  ha  m,  t t re auency  inside  the  COAT 
dither  band,  the  speckle  estimate  t r . > ! . 1 1 . . t same  sine  wave  com- 
bined with  appreciably  reduced  . c ( i utermodulat ion 

products).  With  the  improved  Lty  of  the  estimate 

remained  positive,  indicating  t - .1.1  work  much  better 

than  previously. 

In  our  Initial  exper imen t s , u.  ' i icrformance  of  the  entire 

system  open  and  closed  loop  !■  t > . . me  wave  disturbances: 

at  6.5  kHz  (the  second  lowest  ) and  at  13.4  kHz 

(the  second  highest  COAT  dithei  lie  • results  are  shown  in 

Figures  36  through  41.  Figure  . • t i : i spectrum  of  the 

signal  from  the  target  FMl  me.  < ting  conditions  without 

speckle  interference.  Under  tia  .c  is  well  converged 


The  new  dither  t requeue  is 

15.0  kHz. 


Unless  otherwise  spei  il i 
this  final  configuration 
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Figure  36.  Fourier  spectrum  of  target  PMT  signal 
for  normal  COAT  operating  conditions. 
Horizontal  scale:  2 kHz/division. 
Vertical  scale:  10  dB/division. 


on  Che  target  and  the  spectral  power  at  each  dither  frequency  is  quite 
small.  These  signals  consist  largely  of  noise,  small  amounts  of  energy 
at  the  dither  frequencies,  and  various  intermodulation  produce  fre- 
quencies. When  we  inject  the  single-frequency  sine  wave  plus  dc  offset 
speckle  multiplier  via  the  speckle  simulator  and  at  the  same  time  effec- 
tively open  the  speckle  detector  loop  by  channeling  dc  voltage*  to  the 
denominator  of  the  speckle  divider,  we  disable  the  COAT  channel  having 
that  same  frequency  and  to  some  extent  the  remaining  channels  as  well. 

This  is  evidenced  by  the  increased  power  at  6.5  kHz  (shown  in  Figure  37(a)) 
and  by  the  increased  power  at  13.5  kHz  (shown  in  Figure  38(a)).  In  both 
figures,  there  is  a slight  increase  in  spectral  power  throughout  the  COAT 
dither  band,  which  indicates  lower  overall  convergence.  When  the  ASPECT 
speckle  detection  and  estimation  loop  is  activated,  most  of  the  disturb- 
ing sine  wave  is  cancelled  and  the  power  spectrum  of  the  target  PMT  signal 
is  returned  to  normal,  as  shown  in  Figures  37(b)  and  38(b). 

Figures  39,  40,  and  41  show  the  effects  of  the  above  sine  wave  dis- 
turbances on  the  PMT  signal  as  a function  of  time.  The  mean  value  of 
this  signal  is  directly  proportional  to  target  illumination  and  is 
therefore  indicative  of  the  degree  of  COAT  convergence.  The  top  trace 
in  Figure  39  shows  the  reference  ideal  signal  with  no  speckle  interference. 
The  middle  trace  shows  the  degraded  signal  that  results  when  the  6.5-kHz 
sine  wave  is  introduced  and  there  is  no  ASPECT  cancellation.  The  bottom 
trace  is  just  the  reference  zero  voltage  resulting  from  totally  blocking 
the  target  PMT  aperture.  The  top  trace  in  Figure  40  shows  the  improved 
signal  level  resulting  from  the  ASPECT  cancellation  of  the  6.5-kHz  sine 
wave  modulation.  The  middle  and  bottom  traces  are  the  same  as  those  in 
Figure  39.  Finally,  Figure  41  shows  the  improved  signal  level  and  the 
zero  reference  that  result  from  the  ASPECT  cancellation  of  the  13.4-kHz 
sine  wave.  The  degraded  signal  for  this  case  is  not  shown,  but  it  was 
virtually  indistinguishable  from  the  6.5-kHz  case.  Table  1 summarizes 


This  is  done  by  adjusting  the  corner  frequency  on  the  usual  15-kHz  low- 
pass  filter  to  a much  lower  value  making  the  divider  a simple  gain  con- 
trol in  the  COAT  servo  loop.  Since  there  is  an  automatic  gain  control 
immediately  following  the  divider,  this  modified  gain  does  not  change  the 
normal  COAT  servo  behavior. 


Table  !•  Comparison  of  Experiment  and  Computer  Simulation  Data 


Disturbance 

Normalized  Target 
Glint  Irradiance3 
without  ASPECT 

Normalized  Target 
Glint  Irradiance3 
with  ASPECT 

Sinewave  plus  dc  offset*3 

0.39 

0.91 

6.5  kHz 

Sinewave  plus  dc  offset*3 

0.39 

0.65 

13.4  kHz 

Speckle  noise  of  figure 

0.48 

0.92 

The  irradiance  is  normalized  to  the  no  speckle  noise  level. 
^Ratio  of  ac  to  dc  was  approximately  0.8. 


Fourier  spectrum  of  target  FMT 
with  13.4  kHz  simulated  speckle 
modu 1 at  ion . 
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Target  PMT  signal  versus  time.  Top  trace:  ASPECT 
restored  level  with  6.5  kHz  sinewave  modulation. 
Middle  and  bottom  traces:  same  as  Figure  39. 


Figure  40 
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the  observed  target  illumination  values.  In  the  6.5-kHz  case,  very 
nearly  full  restoration  of  target  illumination  was  achieved  by  the  ASPECT 
speckle  cancellation  system;  in  the  13.4-kHz  case,  only  partial  restora- 
tion was  observed.  This  is  consistent  with  our  computer  simulation 
results  (discussed  in  Section  4.C).  We  also  observed  the  behavior  of 
the  seven-channel  COAT/ASPECT  system  when  subjected  to  the  noise  shown 
in  Figure  35.  With  no  noise  from  the  speckle  simulator,  the  target  PMT 
output  of  Figure  42(a),  which  is  proportional  to  the  target  irradiance, 
was  recorded.  The  denominator  input  to  the  speckle  divider  (i.e.,  the 
estimate  of  the  speckle  modulation  function)  and  the  divider  output  are 
shown  in  Figure  42(b).  Since  there  was  no  noise,  the  estimated  modula- 
tion function  (the  top  trace)  is  a constant.  The  divider  output  (the 
bottom  trace)  is,  therefore,  just  proportional  to  the  target  PMT  signal. 
Figure  43(a)  shows  the  degraded  target  irradiance  resulting  from  the 
introduction  of  the  multiplicative  noise  from  the  speckle  simulator,  the 
ASPECT  speckle  detection  loop  being  inactive.  It  is  0.48  of  the  refer- 
ence no  speckle  noise  level  of  Figure  42(a).  Figure  43(b)  shows  the 
ASPECT  estimate  of  the  speckle  noise  in  the  top  trace  and  the  divider 
output  in  the  bottom  trace.  When  this  esimtate  is  fed  into  the  denominator 
input  of  the  speckle  divider  closing  the  ASPECT  loop,  the  divider  output 
changes  to  that  of  the  bottom  trace.  Figure  44(b).  The  target  irradiance, 
shown  in  Figure  44(a)  is  restored  to  0.92  of  the  no  speckle  noise  level, 
about  a factor  of  two  improvement. 

C.  CONCLUSIONS 

The  experimental  study  reported  here  clearly  shows  that  multiplicative 
noise  can  ba  estimated  and  cancelled  by  a hardware  ASPECT  system.  In 
retrospect,  the  laboratory  hardware  could  have  been  improved  if  provisions 
had  been  made  to  accommodate  a much  higher  ASPECT  probe  carrier  fre- 
quency. But  despite  this  shortcoming,  the  evidence  is  still  quite  con- 
clusive. The  laboratory  hardware  performed  as  well  as  we  had  anticipated, 
based  on  computer  simulation  results  (which  had  indicated  that  reasonable 
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SECTION  6 


AUTOMATIC  BANDWIDTH  ADAPTIVE  CONTROL  (ABAC) 

A.  GENERAL  DESCRIPTION 

Previous  work1’^’1"  has  shown  that  there  is  a direct  correlation 

between  the  amount  of  COAT  performance  degradation  due  to  specile  inter- 

* 

ference  and  the  magni  tude  of  the  system's  servo  bandwidth.  That  is  to 
say,  there  will  be  less  speckle  degradation  at  smaller  bandwidths  and 
vice  versa.  Conversely,  if  one  wants  to  compensate  a dynamic  phase  error, 
such  as  atmospheric  turbulence,  the  opposite  holds  true.  Too  small  a 
bandwidth  will  prevent  the  servo  from  compensating  the  rapidly  changing 
phase  errors.  Clearly,  there  is  some  optimum  value  of  bandwidth,  above  or 
below  which  the  performance  degrades.  However,  it  is  impossible  to  anti- 
cipate what  value  this  should  be  because  it  is  a function  of  the  amplitude 
and  frequency  distribution  of  both  of  the  turbulence  and  the  speckle, 
making  it  different  for  every  scenario.  lhe  ABAC  system  is  designed  to 
actively  and  automatically  seek  this  optimum  bandwidth  for  all  scenarios. 

In  Figure  45,  let  the  signal  S in  m the  COAT  receiver  be  given  by 
Eq.  3.16.  A dither  modulation  of  the  COAT  servo  bandwidth  f at  frequency 
oj  produces  a temporal  modulation  of  the  target  glint  intensity,  given  by 


91 

Ai  = — ~ At  sin(m  t)  , (6.1) 

g <fs  s g 

where  Af^  is  the  amplitude  of  the  bandwidth  dither.  If  there  is  very 
little  speckle  modulation  at  the  dither  frequency,  then,  after  synchro- 
nously detecting  and  low-pass  filtering  signal  S^,  we  get  the  gain  or 
bandwidth  correction  signal,  G(t)  given  approximately  by 


The  term  bandwidth,  as  used  throughout  this  section,  refers  to  the  value 
of  bandwidth  at  full  convergence.  Section  7. A discusses  how  the  band- 
width varies  with  convergence  level. 
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where 

K i optical  detector  constant 

G^  = the  adjustable  gain  parameter  in  the  ABAC  loop. 

Eq.  6.2  illustrates  a fundamental  problem  associated  with  the  ABAC  servo. 
The  derivative 


3f 

s 


cannot  be  anticipated  without  a priori  knowledge  of  the  nature  of  the 

speckle  and  turbulence  degradation  that  will  be  encountered  in  a given 

COAT/target  scenario.  Therefore,  it  is  impossible  to  select  an  optimum 

value  for  the  gain  parameter  G . This  problem  can  be  overcome  by  setting 

g 

the  value  of  G^  very  high  so  that  the  magnitude  of  the  correction  signal 
would  always  be  too  large.  If  a signal  limiter  is  then  placed  after 
the  low-pass  filter  following  the  synchronous  detector,  as  shown  in 
Figure  45,  only  the  polarity  of  the  correction  signal  will  need  to  be 
correct,  because  the  magnitude  will  almost  always  be  the  same.  Thus, 
the  ABAC  servo  will  change  the  gain  in  the  COAT  servo  at  a fixed  rate  in 

either  a positive  or  negative  direction.  The  magnitude  or  clipping  level 
of  the  signal  limiter  should  be  adjusted  so  that  the  gain  does  not  change 
significantly  (e.g.,  ±10%)  during  one  dither  period. 

In  addition,  if  there  is  significant  speckle  modulation  at  the  band- 
width dither  frequency,  then  there  will  be  a spurious  error  signal 
created,  just  as  there  is  in  the  normal  COAT  servo.  There  appears  to  be 
no  way  to  avoid  this  problem  other  than  to  keep  the  ABAC  servo  bandwidth 
low  enough  to  keep  the  speckle  noise  power  below  the  true  signal  power. 
This  would  result  in  a very  slow  servo;  however,  there  are  other  reasons 
for  wanting  a slow  servo,  as  discussed  below. 


121 


It  is  rather  straightforward  to  conclude  that  a very  low  dither 
frequency  for  the  ABAC  servo  (as  compared  to  COAT  dithers)  is  highly 
desirable  if  not  an  absolute  necessity.  Just  for  mechanical  reasons, 
the  dither  period  must  be  longer  than  twice  the  characteristic  response 
time  of  a COAT  servo  (on  the  order  of  a few  milliseconds) . Furthermore 
we  must  recognize  that  the  mean  Strehl  ratio  on  target  from  a speckle- 
or  turbulence-degraded  COAT  system  is  not  well  defined  for  short  periods 
of  time.  It  usually  requires  10  to  20  times  a typical  COAT  convergence 
time  (on  the  order  of  10  msec  or  more)  for  a well-defined  average  conver- 
gence level  to  be  established.  Thus,  for  the  change  in  target  intensity 

3 Ig 

as  a function  of  change  in  bandwidth  (i.e.,  —■£—  in  Eq.  6.2)  to  be  well 

3 f S 

defined,  the  bandwidth  dither  period  must  be  sufficiently  long.  Also, 
because  of  speckle  modulations,  it  will  probably  require  more  time  for 
the  COAT  receiver  intensity  to  establish  a well-defined  mean  level  than 
it  does  at  the  target.  The  conclusion  from  all  of  this  is  that  low  dither 
frequencies  are  most  desirable.  A bandwidth  dither  frequency  of  10  Hz 
for  a 1000-Hz  bandwidth  COAT  servo  is  probably  reasonable. 

Unfortunately,  at  the  beginning  of  this  program  we  did  not  expect 
that  the  dither  would  have  to  be  this  slow.  Only  the  mechanical  consid- 
erations were  anticipated.  Thus,  we  felt  that  a 300-  or  400-Hz  dither 
frequency  and  a 1000-Hz  COAT  servo  would  be  adequate.  The  original  18- 
channel  servo  simulation  was  modified  to  simulate  this  configuration. 

The  nominal  bandwidth  of  the  18-channel  servo  code  ’ is  generally 
between  300  and  500  Hz.  Increasing  the  maximum  allowable  COAT  servo 
bandwidth  to  1000-Hz  would  increase  the  response  time,  thus  allowing 
faster  dithering  of  the  bandwidth,  which  in  turn  would  decrease  the  total 
time  necessary  to  run  each  simulation.  To  implement  this  scheme,  we  had 
to  increase  the  frequency  spacing  between  the  usual  COAT  dither  channels 
from  1400  Hz  to  4000  Hz.  The  18th  dither  frequency  then  became  80.2  kHz, 
which  is  more  than  double  the  previous  value  of  32  kHz.  Consequently, 
to  maintain  the  same  numerical  accuracy  as  before  would  require  us  to 
decrease  the  integration  step  size  by  more  than  a factor  of  2.  Therefore, 
we  chose  to  reduce  the  number  of  channels  from  18  to  12,  thereby  reducing 


the  maximum  dither  frequency  so  that  the  integration  step  size  could  be 
left  approximately  the  same,  6t  = 4 psec.  With  these  modifications, 
combined  with  computer  core  constraints,  we  were  able  to  make  simulation 
runs  lasting  for  192  msec  at  reasonable  costs.  This  allowed  us  to  use 
an  ABAC  dither  frequency  of  50  Hz  and  still  run  simulations  for  almost 
10  dither  cycles.  Even  through  a lower  dither  frequency  and  many  more 
dither  cycles  would  be  much  preferable,  the  data  gathered  under  these 
constraints  is  still  fairly  conclusive. 

A dynamic  phase  error  <f>n(t)  was  also  introduced  into  each  COAT 
channel  to  simulate  atmospheric  turbulence.  This  was  done,  using  the 
following  expression  for  the  phase  error  in  an  arbitrary  channel  n : 


$n(t)  = Atr  si n (wtrt  + Yn),  (6.3) 

where  A , o>  , and  y are  input  parameters.  The  variables  Y_  are  chosen 
tr  tr  n n 

randomly  between  0 and  2ti.  Including  this  simulated  turbulence  gives 
us  a more  realistic  perspective  of  ABAC  operations. 

B.  COMPUTER  SIMULATION  RESULTS 

We  have  run  closed  loop  simulations  of  the  ABAC  servo  system  using 
simulated  turbulence  and  the  same  three  speckle  modulation  functions 
described  earlier  in  Section  4.C.  We  will  refer  to  them  by  specifying 
the  target  rotation  rate.  The  magnitude  of  the  turbulence  in  each  case 
was  A = tt/3  rad,  and  the  frequency  was  1256.64  rad/sec  (200  Hz).  Some 
open  loop  results  are  given  below  to  give  the  reader  a feel  for  what  one 
should  expect  from  an  ABAC  system.  Figures  46,  47,  and  48  show  the 
average  Strehl  ratios  for  the  12-channel  COAT  system  described  in  Sec- 
tion 6. A as  a function  of  servo  bandwidth.  These  averages  were  obtained 
from  computer  simulation  runs  of  the  COAT  system  operating  at  fixed  band- 
width and  subjected  to  the  indicated  speckle  and  turbulence  degradations. 
These  data  show  approximately  where  the  optimum  bandwidth  for  each 
scenario  should  be.  However,  there  is  a great  deal  of  fluctuation  from 
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Figure  47.  COAT  system  performance  at  fixed  bandwidth  targer  rotation  rate 
2.0  rad/sec. 


AVERAGE  STREHL  RATIO 


these  average  Strehl  ratios  at  any  given  bandwidth  for  short  periods  of 
time.  The  error  bars  in  each  figure  indicate  the  standard  deviation. 
These  fluctuations  can  cause  the  average  Strehl  ratio,  over  time  periods 
on  the  order  of  a bandwidth  dither  period,  to  vary  considerably  from  this 
mean  value.  Figures  47  and  48  show  that  the  optimum  bandwidth  for  target 
rotation  rates  2.0  and  10.0  rad/sec,  is  100  Hz  or  less  most  of  the  time. 
But  for  short  periods  of  time,  some  higher  bandwidth  may  be  optimum. 

This  phenomena  is  more  apparent  in  Figure  46,  where  it  is  difficult  to 
determine  .'hat  the  optimum  bandwidth  for  a target  rotation  rate  of 
0.4  rad/sec  really  is.  Since  the  ABAC  servo  responds  on  a time  scale 
much  shorter  than  the  192  msec  over  which  these  data  were  gathered,  these 
fluctuations  do  not  really  present  a great  problem.  The  ABAC  servo 
simply  tries  to  compensate  these  fluctuations  by  changing  the  bandwidth 
in  a direction  that  increases  the  Strehl  ratio.  In  other  words,  ABAC 
does  not  really  seek  a long  time  optimum  bandwidth;  rather,  it  seeks  a 
long  time  optimum  Strehl  ratio,  which  is  really  the  parameter  we  wish  to 
optimize. 

This  can  best  be  appreciated  in  Figures  49  and  50.  Figure  49(a) 
shows  COAT  performance  with  a closed-loop  ABAC  servo  in  the  form  of  tar- 
get Strehl  ratio  versus  time.  Figure  49(b)  shows  the  ABAC  bandwidth 
correction  signal.  As  discussed  above,  only  the  polarity  of  this  signal 
is  important  because  the  ABAC  servo  generates  a constant  change  in  band- 
width in  either  the  positive  or  negative  direction  in  accordance  with 
this  polarity.  The  large  oscillations  at  the  beginning  of  the  plot  are 
system  transients  due  to  the  initial  conditions.  Figure  49(c)  is  a plot 
of  the  average  Strehl  ratio  over  24-msec  intervals  superimposed  on  a plot 
of  ABAC-generated  COAT  servo  bandwidth.  The  bandwidth,  which  was 
initialized  to  600  Hz,  starts  to  decrease  when  the  ABAC  servo  is  turned 
on  40  msec  after  the  start  of  the  simulation.  The  bandwidth  continues 
to  decrease  as  long  as  the  Strehl  ratio  increases  (or  remains  constant). 
When  the  bandwidth  goes  below  300  Hz,  the  Strehl  ratio  decreases  slightly 
and  the  ABAC  system  reverses  the  change  in  bandwidth,  with  a resultant 
further  increase  in  the  Strehl  ratio.  The  final  bandwidth  after  192  msec 
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Figure  49.  COAT  system  performance  with  open  loop  ABAC 
GRC-1.  Initial  bandwidth:  600  Hz. 
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is  almost  700  Hz.  From  Figure  49(b),  however,  it  appears  that  the  ABAC 
servo  correction  signal  was  just  about  to  change  polarity  again  which 
would  have  caused  a reversal  of  the  bandwidth  increase.  It  is  unlikely 
that  the  ABAC  servo  would  have  found  an  optimum  bandwidth  in  this  case. 

It  seems  more  likely  that  it  would  continue  searching  somewhere  between 
200  and  800  Hz.  This  search  would  not  be  random,  however,  because  the 
servo  would  constantly  be  trying  to  maintain  the  highest  Strehl  ratio 
possible  over  any  given  dither  period.  The  open  loop  data  of  Figure  46 
is  quite  consistent  with  this  interpretation  because  there  is  only  a 
small  amount  of  variation  in  average  Strehl  ratio  over  this  range  of 
bandwidth  and  there  is  considerable  fluctuation. 

To  be  sure  that  the  ABAC  system  would  approach  the  same  state  or 
the  same  optimum  bandwidth  from  above  or  below,  we  chose  to  run  each 
scenario  twice,  once  with  the  bandwidth  initialized  at  600  Hz,  and  once 
with  the  bandwidth  initialized  at  100  Hz,  and  then  compare  the  results. 
Figure  50  shows  the  same  data  as  shown  in  Figure  49,  for  a 100-Hz 
initialization.  As  is  the  case  in  Figure  49(c),  Figure  50(c)  does  not 
appear  to  show  a long  time  optimum  bandwidth.  These  two  figures  are  poor 
examples  to  compare  this  way  because  of  the  lack  of  definition  of  the 
final  state;  therefore,  we  proceed  to  cases  where  the  target  rotation 
rates  are  2.0  and  10.0  rad/sec  and  the  comparison  is  more  meaningful. 

These  two  scenarios  are  very  similar  in  their  effect  on  COAT  per- 
formance, and  a judgment  of  ABAC  performance  is  quite  simple  to  make. 

In  Figures  51  through  54  (which  are  formatted  the  same  as  Figures  49  and 

50),  the  ABAC  servo  in  each  case  simply  drives  the  COAT  servo  bandwidth 

* 

to  100  Hz  and  keeps  it  there.  This  result  is  consistent  with  the  open 
loop  data  of  Figures  47  and  48,  where  100  Hz  or  less  is  clearly  the 
optimum  bandwidth.  Starting  the  simulations  from  either  600  Hz  or  100  Hz 
clearly  produces  the  same  final  bandwidth. 


Because  there  Is  always  a plus  and  minus  100  Hz  dither  associated  with  the 
ABAC  system,  100  Hz  is  the  lowest  limit  which  we  allow  the  bandwidth 
to  reach.  Therefore,  in  effect,  the  final  COAT  servo  bandwidth  is  not 
100  Hz  but  varies  between  0 and  200  Hz. 
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Figure  50.  COAT  system  performance  with  open  loop  ABAC 
GRC-1,  Initial  bandwidth:  100  Hz. 
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51.  COAT  system  performance  with  open  loop  ABAC  system.  Case: 
GRC-2.  Initial  bandwidth:  600  Hz. 
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Figure  52.  COAT  system  performance  with  open  loop  ABAC  system.  Case: 
GRC-2.  Initial  bandwidth:  100  Hz. 
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Figure  53.  COAT  system  performance  with  open  loop  ABAC 
GRC-3.  Initial  bandwidth:  600  Hz. 
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Figure  54.  COAT  system  performance  with  open  loop  ABAC 
GRC-3.  Initial  bandwidth:  100  Hz. 
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CONCLUSIONS 


This  section  has  demonstrated  that  an  ABAC  servo  can  optimize  COAT 
servo  bandwidth  and  can  do  this  regardless  of  the  source  of  the  perfor- 
mance degradation.  As  long  as  these  degradations  are  bandwidth-related, 
the  ABAC  servo  will  seek  the  optimum  bandwidth  under  the  given  constraints. 
This  optimization  is  accomplished  by  dithering  the  COAT  bandwidth  and 
then  synchronously  detecting  the  correction  signal  at  the  dither  fre- 
quency. Since  this  signal  is  proportional  to  the  magnitude  of  the  change 
in  the  Strehl  ratio  (or  convergence  level)  in  response  to  a smal  change 
in  bandwidth,  and  since  the  response  of  the  COAT  system  to  such  a change 
in  bandwidth  is  very  slow,  particularly  under  conditions  of  speckle 
interference,  the  bandwidth  dither  must  also  be  very  slow.  We  have 
found  that  the  dither  frequency  should  be  about  two  orders  of  magnitude 
lower  than  the  COAT  servo  bandwidth. 

By  examining  Figures  49  through  54,  one  might  conclude  that  ABAC, 
or  more  generally  bandwidth  variation,  is  only  a marginally  useful  concept 
because  the  improvement  in  Strehl  ratio  in  any  one  scenario  is  not  that 
dramatic.  The  most  impressive  performance  is  shown  in  Figure  53,  where 
the  Strehl  ratio  increases  from  about  0.3  to  0.5.  However,  a realistic 
COAT  target  will  not  be  a sphere  rotating  at  a constant  rate  such  as  the 
ones  modeled  in  these  simulations.  A real  target  would  most  likely  have 
complex  geometry  and  move  about  at  various  and  continually  changing 
velocities  and  angular  velocities.  This  would  produce  a constantly 
changing  speckle  pattern.  This,  combined  with  atmospheric  turbulence 
and  any  other  bandwidth  sensitive  distortions,  would  create  an  ever 
changing  scenario.  It  is  ABAC's  ability  to  adapt  the  bandwidth  to  the 
scenario  in  an  optimum  manner  that  best  illustrates  its  utility. 


SECTION  7 


ANALYTICAL  PREDICTIONS  OF  SPECKLE-RELATED  PERFORMANCE  DEGRADATION 
OF  MULTIDITHER  COAT  SYSTEMS 

This  section  expands  on  an  analysis  originally  developed  under 
contract  D30602-76-C-0021  and  reported  in  Refs.  1 and  12.  In  that 
analysis,  we  were  able  to  calculate  the  time-averaged  convergence  level* 
of  a COAT  system  subjected  to  multiplicative  noise  caused  by  speckle  as 
a function  of  the  relative  noise  energy.  More  precisely,  we  calculated 
the  convergence  level  as  a function  of  a parameter  called  speckle  coeffi- 
cient, or  Cg,  which  is  just  the  square  root  of  the  average  noise  power 
within  the  COAT  dither  frequency  bands.  If  we  assume,  as  we  did  in  the 
original  analysis,  that  the  gain  or  bandwidth  characteristics  of  a COAT 
system  are  independent  of  convergence  level,  then  a knowledge  of  these 
and  the  spectrum  of  the  speckle  modulation  function  is  all  that  is 
necessary  to  compute  a single  value  of  Cg,  which  then  allows  us  to  calcu- 
late the  convergence  level.  In  real  multidither  COAT  systems,  however, 
the  open  loop  dc  gain  varies  with  convergence  level.  Therefore  in  prin- 
ciple there  exists  a separate  value  of  Cg  for  every  value  of  convergence 
level.  This  dependence  of  Cg  on  convergence  level  can  be  thought  of  as 
a constraint  relationship,  whereas  the  original  relationship  which  allows 
us  to  calculate  the  average  convergence  level  as  a function  of  Cg  is  a 
statistical  relationship  based  on  physical  principles.  Both  of  these 
relationships  must  be  satisfied  simultaneously  to  represent  a real  COAT 
system  in  a true  physical  scenario.  Since  the  latter  relationship  is 
already  known  from  Ref.  1,  we  now  proceed  to  derive  the  constraint 
relationship.  To  do  this  requires  calculating  the  open  loop  dc  gain 
associated  with  the  speckle  interference  in  an  arbitrary  channel  n. 

In  Ref.  1,  we  have  shown  that  the  synchronously  detected  signal  in 
each  channel  consists  primarily  of  the  sum  of  two  components:  the  true 


* 

In  this  context,  the  convergence  level  is  just  the  Strehl  ratio  on  a 
target  glint. 
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error  correction  signal  and  the  interference  signal  caused  by  speckle. 
Furthermore,  we  have  demonstrated  that  the  convergence  level  attained 
by  the  COAT  system  was  always  such  that  the  average  power  after  synchronous 
detection  was  equal  in  each  of  these  signals.  Therefore,  we  assumed  that 
the  average  open  loop  dc  gain  associated  with  the  speckle  interference 
is  the  same  as  that  for  the  true  error  signal.  Based  on  this,  we  can 
calculate  the  expectation  value  of  this  parameter  in  an  arbitrary  channel 
n of  an  N-channel  multidither  COAT  system. 

Eq.  11  of  Ref.  1 exprsses  the  normalized  sum  of  all  the  COAT  receiver 
intensity  modulations  at  the  dither  frequencies.  The  modulation  at  the 
n*"*1  frequency  I , is  0ust  one  term  in  that  sum.  Therefore,  1^  can  be  written 


('I'  )Ji  ('I'  ) r q 

r — I / sin  (6  -6  )|  sinto  t, 

2 L z— ' n m J n 

N _ 1 

m=l 

m/n 


(7.1) 


where  the  variable  t is  time,  is  the  amplitude  of  the  dither  modulation 
in  radians  of  wavelength,  Jq  and  are  the  zeroth  and  first  Bessel  func- 
tions, u is  the  dither  frequency,  and  3^  is  the  average  phase  at  the 
target  associated  with  COAT  channel  m.  Using  trigonometric  identities, 

Eq.  7.1  can  be  rewritten  as 


4J.CU-  ) 

i o 


V_  sin  (6  -B_  ) sinui  t, 
Rn  n Rn  n 


(7.2) 


where  V_  and  3„  are  defined  as  the  resultant  amplitude  and  phase  of  the 
Rn  Rn 

f lowing  vector  sums: 


We  will  assume  a COAT  system  without  an  AGO,  as  was  done  in  the 
original  analysis. 


VRn  C0S  (eRn}  = 


V J (ifi  ) cos  (8  ) 
/ , o o m 


(7.3) 


VRn  Sin(BRn) 


= Y J (<l<  ) sin(6  ). 

o o m 


(7. A) 


Physically,  SRn  is  the  average  relative  phase  of  the  resultant  field 

at  the  target  glint  associated  with  all  COAT  cahnnels  except  channel  n, 

and  VRn  is  proportional  to  its  amplitude.  This  can  be  seen  as  follows. 

In  the  original  statistical  analysis  (see  Ref.  1),  we  assumed  that  the 

target  glint  was  illuminated  by  N identical  sources,  each  having  a field 

amplitude  of  N ^ and  a phase  of  6 + U;  sinu)  t,  where  n refers  to  the  n^ 

no  n 

element  or  channel  in  a zonal  multidither  COAT  system.  Since  the  dither 

phase  ip  + sin  w t changes  at  a much  faster  rate  than  the  means  phase  S , 
o n ^ n 

we  may  represent  the  n element  field  at  the  target  by  an  average  phasor 

quantity  V , which  is  the  time  average  of  the  true  field  over  a dither 

period  T : 

n 


- y 

JT  / 

n J 


i($  sinco  t)  J (ij/  ) i$ 

no  n o o n 

dt  e = e 

N 


Except  for  the  normalization  factor  N , the  summations  in  Eqs.  7.3  and  7. A 
are  just  the  component  vector  sums  of  average  fields  defined  this  way. 

Thus,  in  Eq.  7.2  is  just  the  phase  error  in  channel  n that  must  be 

corrected  by  the  COAT  servo. 

The  correction  signal  is  proportional  to  the  amplitude  of  the 
dither  modulation  at  frequency  w . Thus,  using  Eq.  7.2  we  define: 


s 4 


4W 


V„  sin (6  -B  ) 
Rn  n Rn 
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where  is  a proportionality  constant  as  well  as  an  adjustable  gain 
parameter  (see  Figure  13)  and  the  factor  1/2  comes  from  the  synchronous 
detection  process.  The  open  loop  dc  gain  in  channel  n is  defined  as 


We  are  now  in  a position  to  calculate  the  expectation  value  of  g 

n 

based  on  the  statistical  model  developed  in  Ref.  1.  To  do  this,  we  will 
first  derive  the  relationships  among  6^,  vRn>  and  BRn  at  a specified 
convergence  level.  The  total  average  resultant  amplitude  V and  phase 
8 associated  with  the  entire  field  at  the  target  glint,  analogous  to 
the  resultant  field  defined  in  Eqs.  7.3  and  7.4  can  be  defined  as 


The  variable  V can  easily  be  related  to  the  mean  convergence  level  1 
K ] 

defined  in  Eq.  10  of  Ref.  1: 


VR  and  V are  just  the  amplitudes  of  phasors  V and  V, 


This  definition  of  6R  is  inconsistent  with  Ref.  1 
is  identical  to  0Rn  in  this  report. 


(7.12) 


r 


Their  phasor  relationship  is  shown  graphically  in  Figure  55.  It  can  be 
expressed  mathematically  as: 


VR  ' V + W “S<W  <'-l3» 

° ■ V Sln<6R»-6R)  + W sln(B„-V  (7-14) 

Using  Eqs.  7.10,  7.13,  and  7.14,  it  can  readily  be  shown  that 
Eq.  7.7  may  be  rewrittei;  as: 


2J  fo)  G sin(e  -B„)  , , 0 , 1/2 

n = ~ 1 2 2 » \ “ N[1_Jn(^)]  • (7'15> 

n nW  ) 1 M o o ) 


n Rn 


The  phase  angle  B , which  is  the  resultant  absolute  phase  of  the  total 
field  at  the  target  glint,  can  have  any  arbitrary  value.  It  is  conven- 
ient for  computation  to  define  it  to  be  zero.  We  can  then  compute  <g  >, 
the  expectation  value  of  the  open  loop  dc  gain,  consistent  with  the 
statistical  model  of  reference  1,  by  using  the  following  relationship: 


l > = — / d8  g 

bn  a I n n 

a/2 

2 2 1/Z  f 

2W  Gi(N  VN[1-Jo  ^o’11  / 

- J 


sing 


dB 


aN 


-a/ 2 » (en~BRn) 


(7.16) 
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The  factor  a is  the  probability  density  as  defined  in  Ref.  1.  The 

limits  of  integration  in  Eq.  7.16  are  chosen  to  be  consistent  with  the 

choice  of  3 =0.  The  phase  variable  3 can  be  defined  in  terms  of  3 
K Rn  n 

using  Eqs.  7.13  and  7.14  and  the  integral  in  Eq.  7.16  can  be  computed 
numerically  to  give  the  expectation  value  of  g^  at  any  specified 
convergence  level. 

The  COAT  servo  bandwidth  characteristics  are  obtained  by  multiplying 
the  open  loop  dc  gain  by  T(w),  the  transfer  function  of  the  low-pass 
filter  following  the  synchronous  detectors.  Thus,  the  expected  open  loop 
gain  g(w)  is: 


g(w)  = <gn>  T(w). 


(7.17) 


The  average  closed- loop  complex  gain  G(w)  is  then  given  by  the  usual 
relationship: 


G (w)  = 


g(w)  + 1 


(7.18) 


If  desired,  the  average  servo  bandwidth  for  any  specified  convergence 

level  and  C at  that  level  can  be  calculated  as  described  in  Ref.  1. 
s 

A more  accurate  way  of  computing  Cg  is  to  first  calculate  the  average 

power  contributed  by  the  speckle  modulations  in  each  channel  using  the 

complex  gain  profile  G(w)  and  then  to  sum  over  all  channels  to  get  the 

total  average  power.  Thus,  we  present  the  following  modified  definition 

of  C : 
s 


N °° 

1/ 

n=l  °° 


a(u-w  )a  (w-u)  ) G(co)  G (u>) 

n n 


(7.19) 


The  variable  T is  one-half  the  time  period  over  which  the  speckle  modula- 
tions are  def ined  and  over  which  the  mean  convergence  level  is  to  be 
predicted.  The  function  a(w)  is  the  Fourier  transform  of  the  speckle 


143 


modulation  function,  and  is  the  n dither  frequency.  The  frequency 
shift  of  the  function  a(u>)  in  Eq.  7.19  accounts  for  the  synchronous  detec- 
tion process.  Even  though  the  limits  of  integration  in  Eq.  7.19  are  ±°°, 
the  values  of  G(w)  beyond  several  servo  bandwidths  are  typically  so  much 
lower  than  those  near  zero  frequency  that  an  adequate  numerical  value  of 
Cg  may  be  obtained  by  integrating  over  plus  and  minus  several  bandwidths. 

A computer  program  which  calculates  Cg  as  a function  of  convergence 
level  in  the  manner  described  above  has  been  coded.  This  program  also 
calculates  the  predicted  convergence  level  as  a function  of  Cg  based  on 
the  original  analysis  of  Ref.  1.  A typical  program  output  is  shown  in 
Figure  56.  The  curve  marked  physical  relationship  corresponds  to  the 
analysis  of  Ref.  1,  while  the  one  marked  constraint  relationship  corre- 
sponds to  the  present  analysis.  The  mean  Strehl  ratio  or  convergence 
level  where  the  two  relationships  have  simultaneous  solution  (i.e.,  where 
the  two  curves  inversect)  is  the  analytically  predicted  value,  about  0.45 
in  this  case. 

In  both  of  these  computations,  we  have  assumed  a 12-channel  multi- 
dither COAT  system  having  a dither  amplitude  of  = 0.349  rad  of  wave- 
length in  each  channel.  The  dither  frequencies  range  from  8.0  to  52.0 
kHz  separated  by  intervals  of  4.0  kHz.  In  calculating  Cg,  we  used  the 
speckle  modulation  function  from  Figure  12  in  Ref.  1.  The  transfer 
function  T(w)  was 

-1  -4 

T(w)  = (1  + i— ) (1  + i-f)  , (7.20) 

u>0  “i 

where  = 20  irrad/sec,  and  = 1000  irrad/sec  (representing  the  five- 
pole  low-pass  filter  which  follows  the  synchronous  detector  in  each  COAT 
channel).  The  loop  gain  parameter  was  set  at  14,830.  Using  Eq.  7.17, 
this  corresponds  to  an  open  loop,  unity  gain  servo  bandwidth  of  400  Hz 
at  full  convergence.  We  refer  to  this  as  the  nominal  servo  bandwidth. 


r- 
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Figure  56.  Typical  analytical  output.  (This  case  predicts  degrada- 
tion of  lines  from  the  speckle  modulations  of  Figure  12, 
Ref.  4.) 
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These  calculations  were  repeated  at  six  different  values  of  nominal 
servo  bandwidth  for  each  of  the  three  speckle  modulation  functions  of 
Ref.  1,  which  are  the  same  three  functions  that  we  have  been  using 
throughout  this  study.  Figures  57,  58,  and  59  show  the  results  of  these 
calculations.  These  predicted  values  of  convergence  level  are  compared 
to  computer  simulation  results  shown  in  the  same  figures,  representing 
the  same  COAT  servo  configurations  and  degraded  by  the  same  speckle 
modulation  functions.  The  computer  simulation  results  are  time  averaged 
target  Strehl  ratios  or  convergence  levels  taken  over  a period  of  8 msec. 
The  comparison  between  these  results  are  good,  especially  considering  the 
complexities  of  the  real  situation.  The  gain  in  each  channel  varies 
instantaneously  with  the  phase  error  in  that  channel  and  with  the  state 
of  convergence  of  the  remaining  channels.  Thus,  the  channels  always 
have  different  and  varying  bandwidths.  Consequently,  varying  amounts  of 
speckle  interference  are  allowed  to  pass  through  each  channel. 

Also,  once  a channel  has  been  completely  disabled  by  speckle  noise 
interference,  further  increases  in  noise  power  resulting  from  larger 
bandwidths  will  no  longer  cause  more  degradation  of  performance  of  that 
channel.  Furthermore,  if  the  spectral  power  of  the  speckle  modulation 
function  is  concentrated  in  the  low-frequency  channels,  as  is  the  case 
in  Figures  57  and  58  (and  to  some  degree  in  all  speckle  modulation  func- 
tions), then  increasing  the  bandwidth  will  not  cause  an  equally  distri- 
buted increase  in  the  speckle  interference  among  all  channels.  Most  of 
the  increase  will  occur  in  the  low-frequency  channels.  After  these  low- 
frequency  channels  are  completely  disabled,  then  further  increases  in 
bandwidth  will  cause  a proportionately  lower  decrease  in  overall  COAT 
performance.  This  is  probably  the  reason  for  the  consistently  greater 
discrepancy  between  analysis  and  computer  simulations  results  at  the 
higher  nominal  bandwidths.  The  computer  simulations  reflect  the  true 
situation,  while  the  statistical  analysis  which  treats  all  channels  as 
having  the  same  average  phase  error  related  to  the  total  average  power 
cannot  account  for  this  phenomenon. 
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Figure  57.  Comparison  of  analysis  and  computer  simulations  using  speckle 
modulations  from  Figure  10,  Ref.  1. 
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Figure  58.  Comparison  of  analysis  and  computer  simulations  using  speckle 
modulations  from  Figure  11,  Ref.  1. 
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Figure  59.  Comparison  of  analysis  and  computer  simulations  using  , 

speckle  modulations  from  Figure  12,  Ref.  1. 
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APPENDIX 


SYSTEM  SAFETY  ANALYSIS 

The  contract  statement  of  work  requires  that  a system  safei  y 
analysis  be  performed  in  accordance  with  paragraph  5.8  2.1  of  MIL -STD-  882. 
Such  an  analysis  was  performed  during  the  contract  work  to  identity  did 
where  necessary,  correct  possible  hazards  associated  wi i h the  >pei at l n 
of  the  system,  known  as  the  "multidither  COAT  experimental  model"  Ihere 
after  called  "the  COAT  system"),  used  in  the  experimental  casks  of  o : i 

k 

contract.  The  analysis  performed  follows  a similar  analysis  pe. 
for  the  same  COAT  system  when  it  was  located  in  a different  area. 

There  are  two  potential  hazard  sources  for  users  cf  this  system 
within  the  laboratory  environment,  where  this  contract  work  was  peifoitned: 
(1)  electrical  shock  and  (2)  eye  damage  by  laser  radiation.  Appropriate 
safety  precautions  were  taken  as  part  of  the  Hughes  Industrial  Safely 
Program  to  ensure  sufficient  fire  protection,  clearance  for  egress  t,on> 
the  laboratory,  access  to  power-line  breakers,  etc.,  so  that  the  require- 
ments of  the  California  OSHA  codes  were  satisfied. 

Electrical  shock  from  the  COAT  system  was  assessed  to  be  in  Cate 
gory  I — Negligible  as  defined  in  paragraph  3.14  of  MIL-STD-882  Proper 
design  of  the  electrical  parts  of  the  system  ensures  that  no  personnel 
can  encounter  lethal  or  hazardous  electrical  voltage  and/or  current 
levels. 

The  hazard  level  for  optical  eye  damage  was  assessed  to  be  in  Cate- 
gory II  — Marginal.  The  nature  of  the  program  required  working  with  tw-.. 
laser  sources,  a 30-mW  He-Ne  unit  operating  at  0.6328  urn  and  a 1 5 W 
argon  unit  with  an  output  wavelength  of  0.488  pm.  Various  system  opt  Us 
reduce  the  maximum  potential  power  levels  to  10  mW  at  0.6328  |im  and  tOO  <nW 
at  0.488  pm.  The  potential  hazard  is  adequately  controlled  by  providing 


* 
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safety  glasses  for  all  personnel  while  working  on  or  around  the  system 

and  by  using  masks,  optical  stops,  etc.  to  confine  the  optical  laser 

beams  to  areas  where  no  personnel  can  directly  view  them  at  day  point 

2 

where  the  optical  power  density  exceeds  2.5  mW/cm  . 

Sufficient  safety  procedures  have  thus  been  implemented  to  ensure 
fhat  both  the  0.6328  pm  radiation  (a  Class  III  hazard  at  the  laser  output) 
and  the  0.488  pm  radiation  (a  Class  IV  hazard  at  the  laser  output)  are 
adequately  controlled  to  reduce  the  hazard  to  Category  II  for  personnej 
working  within  the  laboratory  where  the  COAT  system  is  operated.  In 
addition,  personnel  are  instructed  in  proper  safety  procedures  while 
within  the  laboratory.  Appropriate  warning  signs  notify  visitors  to  tb-. 
laboratory,  and  laser  safety  goggles  are  available  for  visitors. 
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